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Liquid state NMR quantum information processors are dominated by the problem 
of preparing pure initial states. Traditionally these devices employ pseudopure states 
instead of truly pure states, but this approach is not scalable. Furthermore, almost all 
previous demonstrations of information processing in liquid state NMR used separable 
states which can in principle be described by classical models. This thesis addresses the 
two-fold problem of initialization and separability by preparing pure quantum states 
lying above the entanglement threshold. Our pure state quantum computer derives its 
purity from the highly polarized nuclear spin states in the para-hydrogen molecule. 

The thesis begins with a critique of conventional NMR based quantum information 
processing in Chapter Q outlining the major strengths and weaknesses of the technol- 
ogy. Chapter El describes the enhanced magnetic ordering of the nuclear spin states in 
para-hydrogen and Chapter El describes an initialization experiment exploiting this ef- 
fect to achieve pure, entangled states. These states can indeed be used as initial states 
in implementing quantum algorithms. This is shown in Chapter 01 which describes an 
implementation of the Deutsch and the Grover quantum algorithms. 

The "twirl" operation converts a completely arbitrary input state to a Werner 
singlet. The NMR implementation of this operation is taken up in Chapter Chapter 
IH1 analyzes the possibility of sharing the purity of some highly polarized qubits in 
a quantum computer onto quantum subspaces of arbitrary dimensions, and whether 
these sharing operations increase or decrease the likelihood of entanglement. 

Finally, this dissertation has six Appendices which contain notes on singlet-triplet 
mixtures, isotropic mixing sequences, systematic errors in singlet detection, data pro- 
cessing and decoherence modelling. 
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Chapter 1 

Liquid state NMR quantum 
information processing — a critique 

Liquid state nuclear magnetic resonance (NMR) [TJ |21 EH IH E| is a vital spec- 
troscopic technique in analytical chemistry. It also plays a crucial role in present day 
medicine as a valuable imaging technology. Recently, NMR has also sought new rela- 
tionships: interfacing with myriad fields in the biological and physical sciences, probing 
living and non-living specimens of all kinds, and now, in the form of functional imag- 
ing, it is taking the first few strides towards an understanding of the way humans 
think and act! One such exciting partnership is with quantum information processing 
(QIP) jHlini, now blossoming into an established discipline in its own right. This thesis 
explores some facets of this 'NMR-QIP' interface. 

NMR-QIP capitalizes on fifty years of experience in probing nuclear spins in liquid 
solution with radio-frequency (RF) pulses. In fact, several essential features of exper- 
imental QIP, such as Rabi oscillations JU] i n single qubits, coherent manipulations 
of multi-qubit systems and decoherence control have long been developed and imple- 
mented in traditional NMR, decades before QIP was investigated as an independent 
discipline. Today, liquid state NMR is not only a leading technology for demonstrating 
quantum computing (QC) in small physical systems, but also enriches QIP, by inspiring 
new ideas for theoretical and experimental investigation. 

Like all the different ideas being pursued for the physical realization of quantum 
computers, liquid state NMR comes with its strengths and weaknesses. A useful set of 
criteria was laid out by DiVincenzo ^1] and since then, it has become customary to 



2 



evaluate practical QC schemes O HS1 EE! on the basis of these rules, which are listed 
below: 

1. there must exist well-defined, physical qubits and the system should be scalable 
to higher dimensions; 

2. the qubits must be initializable into a well-defined quantum state; 

3. the decoherence times must be long compared to the times for gate operation; 

4. there must also exist a set of quantum gates capable of implementing universal 
quantum dynamics; and 

5. it should be possible to perform single-qubit measurements. 

Most of my D.Phil, research revolves around one of the major problems with liquid 

state NMR 1 , namely initialization (criterion 2). Consequently, in this chapter, I choose 

to present a detailed description of the intrinsic weaknesses in conventional NMR, and 

only touch upon the strengths. Furthermore, this chapter is not intended to become 

another review of NMR-QIP; there are already excellent tutorials and monographs on 

the subject, such as |Hl H3 EH HZl UH1 EH] • Instead, the specific aim of this chapter is to 

introduce the basics, to guide the reader through the material in subsequent chapters 

and also, to come up with a motivation for the forthcoming chapters. 

Liquid state NMR quantum computers, reaching seven qubits [211] , are known to be 

far ahead of competing technologies, at least in terms of the implementation of quantum 

algorithms |20l I2H 1221 123 1231 123 I2H1 - (Recently, the Deutsch-Jozsa algorithm 

has been implemented on a two qubit ion trap quantum computer jHII|.) This relative 

success has been made possible by the ease as well as robustness with which unitary 

operations can be implemented in NMR. These strengths are mentioned in Section fTTTl 

which describes the NMR states, coherences and universal unitary dynamics (criteria 

1, 3 and 4) and also introduces the necessary notation. The problems of initialization 

and single qubit read-out (criteria 2 and 5) are discussed in Sections 11.21 and 11.31 

Another problem, not directly encompassed in DiVincenzo's criteria, is the issue of 

separability and entanglement. Traditional NMR is critically analyzed in view of this 

requirement in Section 11.41 and finally, general methods for addressing the problem of 

low polarizations are addressed in Section 11.51 

1 From here onwards, NMR means 'liquid state NMR', unless the solid state is explicitly mentioned. 
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1.1 NMR States, Coherence and Dynamics 



A qubit is a two level quantum system. A spin s nucleus has 2s + 1 energy levels; the 
qubits in NMR are (generally) spin 1/2 nuclei, and so I often use the words 'qubit' and 
'spin' interchangeably. The two energy levels are labelled |a) and \(3) (or |0) and |1)), 
and these states also serve as the computational eigenstates ^3]. The most readily 
available two level systems are 1 H, 19 F, 31 P, 13 C and 15 N, and what characterizes 
these qubits are their Larmor frequencies |3] , which in the presence of a magnetic field, 
B , are, 

w = ~lB (1.1) 

7 being the gyromagnetic ratio of the nucleus. The qubits are therefore "identified" 
in frequency space, and can be perturbed by applying resonant electromagnetic radi- 
ation. For the nuclear spin levels, the Larmor frequencies (and therefore the resonant 
radiation) are in the RF region of the spectrum. An NMR spectrometer is a highly 
sophisticated, computer controlled radio receiver and transmitter, capable of delivering 
electromagnetic radiation to nuclei in a liquid solution placed in a strong superconduct- 
ing magnet, and of detecting a time-resolved RF signal generated by the nuclei. The 
signal is Fourier transformed and a frequency spectrum is obtained, with the qubits 
identifying themselves by their individual frequencies. In short NMR furnishes well 
characterized qubits, partially satisfying the first DiVincenzo criterion. About the 
other half, that is scalability, I will say more in Section 11.21 

1.1.1 Ensemble quantum computing 

In the NMR setting, each qubit resides in a molecule, along with other nuclei, with 
arbitrary spins. In the sample, there will be a large number of these molecules; for 
example, in 1 ml of a 1 mM solution, there will be roughly 10 17 , rapidly tumbling in the 
liquid. As there is one n qubit quantum computer per molecule, so the sample will con- 
tain an ensemble of quantum computers. All these quantum computers are identical, 
act in parallel and collectively participate in producing a time-dependent macroscopic 
magnetization, which is manipulated and detected in the NMR experiment. A con- 
venient language to describe such ensembles is the density operator formalism [51 13 Ij . 
An n qubit state, for example, is represented as an 2 n x 2 n density matrix, expanded 
in some operator basis. In this thesis, I shall use the notation N = 2 n for the Hilbert 
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dimensions of an n qubit quantum system. Different basis sets, each comprising 4 n 
basis elements, are in use pj; however, one set which is the most common in the NMR 
community is the product operator E2 HUH Ell se ^ which I now describe. 

1.1.2 Product operator basis 

The product operator basis elements are represented as N x N matrices and are sim- 
ply re-scaled versions of products of the Pauli operators [Sj. For a single qubit, the 
correspondence 2 is, 



I2 _ / 1/2 

2 I 1/2 



ko, (1.2) 



\o z . (1.5) 



1 1 /A 

h= = \a x% (1.3) 

1/2 0) 

! ~V 2 \ , 

4 = = hcr y , and (1.4) 

\i/2 y - 

1/2 
-1/2 

The symbol 1 at is an example of a general notation I shall use for anJVxiV identity 
matrix. The matrix, I2/2, is commonly denoted as \E in the NMR literature, but 
I shall keep to the notation introduced above. A single qubit density matrix, p can 
always be uniquely expanded in terms of these four basis elements, 

P = y + E c ^ ( L6 ) 

i=x,y,z 

which has unit trace matrix as it should [S], and the coefficients are given by, 

d = Tr{ph)/Tr{lf) = 27Y(pI i ), (1.7) 

Tr(. . .) representing the trace operation. It is traditional to represent operators of 

the qubits with labels /, S, R, and so on 3 . For a two qubit system, the 4 2 = 16 

2 An operator basis is representation-independent whereas "writing" the basis elements in matrix 

form requires a representation basis. All matrices in this thesis are written in the so-called Zeeman 

basis which is formed by taking the direct products of the Zeeman eigenstates, which are also 

the computational eigenstates. The ordering that I follow in writing these matrices is the traditional 

\a(3), \(3a, |/3/3))} order. 
3 Sometimes, the notation ii, I2 ■ ■ ■ is used, especially for homonuclear systems. 
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product operators are constructed by taking the tensor products of single spin operators 
according to, 

2x Ii<g>Sj, i,j — x,y,z (1.8) 
1 2 

2x — <g> Sj, j = x,y,z (1.9) 
1 2 

2x Ii (g) — , i — x,y,z and (1-10) 

2x 1^®^ = ^. (1.11) 
2 2 2 v ; 

It is also common practice to drop the multiplication signs, and by slight abuse of 

notation, represent the two qubit operators with the same symbols as used for the 

single qubit (e.g., 2 x l 2 /2 ® Sj = 1 2 <8> Sj = Sj), the exact meaning being apparent 

from the context. For example, the sixteen product operators for a two qubit system 

are given in Table ITU Any two qubit state can be expanded in terms of these product 

operators, 



2 [ 2 ' 



i>= o(— )+ CiIi+ Yl diSi + Yl e v 2I * s j- ( Ll2 ) 



i=x,y,z i=x,y,z t=x,y,z 

The factor of 1/2 in front of the scalar element I4/2 ensures that p has unit trace and 
the coefficients c, d and e are calculated in a manner similar to (|1.7|) . Furthermore, the 
I4/4 term does not evolve under unitary operations j3T] and from the point of view 
of the state dynamics, is totally "dormant" . The state is called the maximally mixed 
state and is often dropped out from the state expansion, ()1.12|) . As a result, only 
the traceless operators are written in the expansion and a unit trace scalar operator 
is assumed. The traceless part is called the deviation matrix [7J EI]. Although the 
distinction between the "real" and the deviation state is immaterial for traditional 
NMR and signal detection, for QIP calculations it is better to carefully appreciate this 
difference: therefore, I strictly use the symbol '~' whenever I write only the deviation 
part of the state expansion. Moreover, I represent diagonal matrices with the notation 
{{...}}, with the diagonal elements inside the braces and reserve single braces {. . .} for 
writing general lists of numbers or vectors. For the physical meaning (or lack thereof,) 
of the product operators, one may refer to [3J. 

1.1.3 Coherences 

In Table fTTTl I have grouped together the product operators in terms of their coherence 
order. The concept can be defined in terms of the positions of non-zero elements in their 
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Scalar element I4/2 
Single quantum coherences 



T T 9 <7 

J D "-^a;; 

2I X S Z , 2I y S z , 2I Z S X , 2I z S y 
Izi S z , 2I Z S Z 



Populations 

Multiple quantum coherences 2I X S X , 2I y S y , 2I x S y , 2I y S x 

Table 1.1: The 16 product operator basis elements for a two qubit (IS) system. The 
basis elements are grouped according to their coherence orders. 

matrix representations. The single quantum coherences are located in the positions, 

/ □ □ \ 

□ □ 

□ □ 
□ □ 



1.13) 



/ 



and the operators {I x , I y , S x , S y , 2I X S Z , 2I y S z , 2I Z S X , 2I z S y } comprise only these coher- 
ences. The zero and double quantum coherences are in the positions, 

/ \ / D \ 



□ 



□ 



V 



and 



/ 



□ 



;i.i4) 



/ 



and are collectively designated multiple quantum coherences. Similarly, the diagonal 
terms are called the populations and in certain respects behave like zero quantum terms. 
(The matrix representations of the 16 product operators are tabulated in [SJ |7j.) 

The spin operators involving multiple coherence terms are also expressed in terms 
of their linear combinations, 

(1.15) 
(1.16) 
(1.17) 
(1.18) 



ZQ X — 


2 {2I X S X + 2I y Syj , 


ZQy = 


~^(2I y S x 2I x Syj, 


DQ X = 


- (2I X S X - 2I y S y ) and 


DQy = 


— (2I y S X + 2I x Syj . 
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This allows us to define a new basis, which I call the product operator-multiple quantum 
basis, spanned by the product operators for single quantum terms and populations and 
{ZQ X , ZQ y , DQ X , DQ y } for the multiple quantum terms. 

Two other bases will be occasionally used in parts of this thesis. One is the spherical 
basis |2], in which the single qubit operators are related to the product operators 
through, 

r k =i x ±u y , (1.19) 

and I z and I2/2 are left unchanged; the multiple qubit operators are formed analogously 
to (jl.8|) - (jl.ll|) . by taking tensor products. The other basis I shall occasionally refer 
to, is the polarization operator [2] basis, in which I x and I y are unchanged, and I z and 
I2/2 are replaced by, 

la = Y + h alld (L20) 

h = y' 1 " (L21) 



which are in fact equal to the projection operators |a)(o;| and \(3){f3\. 
1.1.4 Detection in NMR 

All NMR experiments can be described in terms of transfers of coherence orders [HHH]. 
Even the detection process can be clarified using the language of coherences. For 
example, in a standard NMR experiment, only the single quantum coherence terms 
can be detected. As such, they are called the direct observables. Suppose we wish 
to detect a state p: then only the single quantum terms in the expansion (|1.12|) will 
contribute to the signal, while all other terms will remain unobservable, 

P > \Ixi Iy> S X 1 Sy, 2I X S Z , 2IyS Z , 2I Z S X , 2I z Sy\. (1.22) 

The terms {I x , I y , S x , S y } represent in-phase and {2I x S z ,2I y S z ,2I z S x ,2I z S y } indicate 
anti-phase magnetization. The single quantum coherences are also assigned |3] a pos- 
itive or negative sign, depending on whether they are positioned above or below the 
diagonal. The way quadrature detection [3J ESj is conceptually employed in modern 
spectrometers, measures only one kind of, (say, the negative) coherence terms; the 
mechanism is carefully spelled out in 
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The spherical and polarization operators are used in predicting the NMR spectrum 
acquired from a state p. I define a signal vector, Sg, with the terms, 



The two leftmost entries indicate resonances of the I spin, and the two rightmost 
entries indicate resonances of the S spin. Mathematically, for an exponentially decaying 
sinusoidal time domain signal, called the free induction decay (FID), these resonances 
appear as Lorentzian lineshapes [3J EI m the spectrum; real values in Sg indicate 
absorptive and imaginary values represent dispersive character. 

1.1.5 Chemical shift 

The / and S spin resonances, appearing in the NMR spectrum are distinct because the 
/ and S qubits are different nuclei having different Larmor frequencies, (jl.lj) . This is 
most obviously true for a heteronuclear system, involving spins with different chemical 
identities, such as X H and 13 C. However, NMR spectra reveal that even if the / and 
S spins have the same chemical identity, i.e., they form a homonuclear system, their 
frequencies can be different. Each spin is surrounded by a cloud of electrons, which 
shields the externally applied field, B , changing the effective field 4 seen by the nucleus 
to Bq{1 — a). The Larmor frequency, (jl.lj) . of each spin is, therefore, shifted, 



If different spins belonging to a homonuclear system, have different electronic envi- 
ronments and therefore, different chemical shifts, they can still act as distinct qubits. 
A useful parameter is the difference in frequencies, 2tt5 (5 is assumed to be in Hz), 
signifying how far apart the qubit resonances are: for a pair of homonuclear spins such 
as two protons, 5 is of the order of a kHz, whereas for heteronuclear systems, it can 
easily be a few hundred MHz. 

1.1.6 Spin-spin coupling 

For purposes of QIP, it is important, that we have two or more qubits per quantum 
computer, and that these qubits interact with one another. The interaction between 
4 Somc authors use the 'deshielding convention', in which B is replaced by B (l + S). 



{Tr(pl + S a ), Tr(pl+Sp), Tr(pI a S + ), Tr(pIpS+)}. 
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spins can assume different forms [21 HSU EZj , isotropic and anisotropic, direct and in- 
direct, depending on the molecular motions, and the physical state of the sample. In 
solids and liquid crystals, the interactions are dominated by direct, dipole-dipole mag- 
netic interactions between the nuclei. The NMR of such anisotropic media is a separate 
realm altogether; a very good textbook on the subject is and examples of QIP with 
dipolar coupled qubits are presented in |3E1 HUH HO] • 

In liquids, with which we are concerned, motional averaging (Sj leaves us with a 
non-zero isotropic value of a kind of indirect coupling between the spins, which is called 
the J or scalar coupling. This magnetic spin-spin coupling is mediated by electrons 
bonding the nuclei, and is a kind of Fermi-contact interaction. Typical values of J are 
a few Hz for 1 H- 1 H couplings and slightly larger, that is between 100-200 Hz for one 
bond 13 C- 1 H couplings. The signal vector, (J1.23|) suggests the excitation of four lines 
in the spectrum of a two qubit system: in fact, the scalar coupling causes each peak to 
split into two halves. The first two elements, Tr(p I+S a ) and Tr(p I + Sp), for example, 
represent the two / spin resonances, when the neighbouring S spin is in the |a)(a| or 
\P)(P\ state. 

In the NMR tradition, a pair of spins is classified as weakly coupled, if 5 ^> J; the 
system is said to be of the type AX. This will be the case for heteronuclear systems 
as well as for many homonuclear systems in the presence of strong magnetic fields, 
Bq. The system is classified as strongly coupled |7j, and of the type AB, if 8 < J. 
Such systems are the prevalent case in weak or near zero fields. The extreme case is 
an A 2 system, in which 5 = 0, where the spins become equivalent and have identical 
frequencies. Even two heteronuclear spins, in zero field conditions, would form an A 2 
system. 



1.1.7 The NMR Hamiltonian 



For a two qubit, AB system, the NMR Hamiltonian PA^EUE] is given by, 



H = h(u!l z + uj s S z + n J 2 1 • S) 



(1.25) 



where I • S is short-hand for 



(1.26) 
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and in the weak coupling limit, i.e., for an AX system, it becomes, 

H(uj i I z + uj s S z + 7tJ2I z S z ). (1.27) 

I call these Hamiltonians AB and AX Hamiltonians. In these expressions, luj and uig 
represent the chemically shifted Larmor frequencies, (jl.24|) . of the two qubits, and I 
subsequently follow the convention of setting h = 1 . It is also customary, to write the 
rotating frame [3] versions of these Hamiltonians, replacing each u by Q = u — u r f, 
u r f being the frequency of the rotating frame 5 , generally taken to be equal to the 
transmitter frequency |41j . In this dissertation, I shall mostly deal with homonuclear 
systems in the weak coupling limit. 

The Hamiltonian discussed so far is the background or the internal Hamiltonian, 
under which the system evolves at all times. A resonant magnetic field is applied to 
perturb the spins. This control field, B\ is perpendicular to the static field, Bq and oscil- 
lates in the perpendicular plane (the transverse plane) with a frequency ujrf, expected 
to be close to the Larmor frequencies of the spins in order to satisfy the resonance 
condition. The oscillatory field, in addition to a frequency and amplitude parameter, 
also possesses a phase factor <f> in the transverse plane; the control Hamiltonian for a 
perturbation on a single qubit is given by, 

H RF = Ui (l x cos (u RF t + 4>) + I y sin {u RF t + </>)) , (1.28) 

where uo\ = — jBi/2 is called the nutation frequency. The nutation frequency (ampli- 
tude), phase, radio frequency and the time for which this control field is applied, can 
all be modulated in time to achieve the desired state dynamics. 

1.1.8 Dynamics 

Corresponding to a Hamiltonian that is piecewise constant in time, Hi, H2, . . . Hk, 
acting for time intervals ti, t 2 , ■ ■ ■ tk, we can construct a propagator, 

U = T k exp(-iH k t k ), (1.29) 

where T is the Dyson time-ordering operator and a state p evolves according to 
the Liouville-von Neumann equation j2], 

p -> UpU ] ; (1.30) 
5 The symbol for rotating frame £ rf" should not be confused with 'RF' for radio frequency. 
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the correspondence between this equation and the Schrodinger's equation is analogous 
to the relationship between the the density matrix and the pure state. What is referred 
to as a pulse in this thesis, is the pulse propagator, 

P = exp{-iH RF t p ). (1.31) 

Selecting suitable times t p , it is possible to achieve desired flip angles; for example, if 
LUitp = 7r/2, the pulse flips the magnetization through 90°. I shall represent a pulse 
with a flip angle of 0°, with a phase angle by the symbol 0^, and for RF fields along 
the cardinal axes, I shall simply represent the phases by {x,y, —x, —y}- For example, 
in my notation, a pulse with flip angle 180° and phase of 180°, would be written as 
180_2;. So far, I have assumed all axes of rotation to lie in the transverse plane. This 
is not a strict assumption, as rotations around arbitrary axes in the Bloch sphere jS] 
can also be implemented using composite rotations jU EJ H21 HH HH EE] • This includes 
rotations around the ±z axes 

The product operators also facilitate predicting the dynamics of the spin operators. 
We make use of a theorem [3IZj: for a set of three operators satisfying [A,B] = iC 
and its cyclic variations, A evolves under the propagator exp (— z0C) as, 

A exp { ~ i0C \ A cos (0) + B sin (0) . (1.32) 
1.1.9 Quantum logic gates 

The standard network model of quantum computation jH] assumes a network built up 
of logic gates |17j. In NMR, these gates are generally implemented |38| 119*1 loTH I5T] lo^] 
through a sequence of pulses and intervals of free precession (under the internal Hamil- 
tonian). Stated simply, single qubit gates can be implemented by pulses exciting only 
one qubit in the system, whereas two qubit gates normally exploit the J coupling be- 
tween the spins, and often make use of the spin-echo pulse sequence jHHHIHI- Logic gate 
design often involves "sculpting" the Hamiltonian, making use of the concepts learned 
from average Hamiltonian theory j2|. Sometimes, use is made of frequency selective 
pulses, so called soft pulses exciting peaks from a single qubit, or "finer tuned", 
transition selective pulses exciting only one peak in a group. These selective 

pulses, also come with their concomitant problems of cross-talk with other qubits and 
coupling evolutions during finite pulse lengths; different methods have been devised 
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EI3 EZ] to ameliorate these effects. In our work [5B1 EH] , we mostly implemented 
selective pulses using hard pulses and periods of free precession, a technique based 
on the Jump-and- Return jHO] sequences. I shall return to this method in Chapter El 
Furthermore, in Chapter we achieve selective excitation, using a method devised by 
Fortunato et al. which numerically searches for a shaped pulse profile, that max- 
imizes some pre-defined fidelity measure. Experimentally, all pulses, be they soft or 
frequency non-selective {hard) 4J, always have instrumental imperfections leading to 
systematic errors [TBI EH] , and these can also be tackled [U21 EH EH EH] very effectively. 
(These errors are also considered in Appendix ??.) In this dissertation, I represent soft 
pulses achieving a flip angle of 9, as 01^ or OS^, exciting the / or S spin, and hard 
pulses as 9$ without a spin label. 

It was shown early on j66j that a controlled-NOT (c-NOT) gate together with single 
qubit rotation, suffices for implementing universal quantum computation jHZI- NMR 
readily provides c-NOT and single qubit rotations, and can therefore, simulate univer- 
sal dynamics, to some required accuracy, satisfying the fourth DiVincenzo criterion. In 
fact, a major factor contributing to the early success of liquid state NMR, was the reli- 
ability and flexibility of applying arbitrary unitary propagators to the system, making 
it a convenient test-bed for many quantum experiments. 

1.1.10 Decoherence 

Decoherence jS] is the decay of coherences, leading to the gradual decay of off-diagonal 
terms in the density matrix. The time-scale over which this decay takes place is im- 
portant, because quantum computations can only last for times shorter than this char- 
acteristic length of time, putting an upper limit on the length of computations, before 
the quantum computer "collapses" into a set of so-called pointer states [HE] • In NMR, 
decoherence is captured by the term relaxation PjJ El El E] and two characteristic time 
scales are defined: T\ which determines the return of the state to the equilibrium con- 
dition and T 2 , which accounts for phase randomization 6 . (A detailed treatise on the 
theory of relaxation in NMR is [69|.) For all spin 1/2 nuclei, T% < T\, and therefore, 
in a simple decoherence model, T 2 is assumed to be the limiting coherence time scale. 

Typically, a two qubit gate can be completed in a time proportional to 1/ J s, so a 
6 Strictly speaking, the T\ relaxation also contributes to the phase randomization. 
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rough measure of the number of (two qubit) operations that can be completed in the 
coherence time is T2/(l/J), a parameter we desire to maximize. I shall return to this 
argument in the introduction to Chapter 0] It is difficult, in general, to estimate the 
number of two qubit operations that can be performed before decoherence takes its 
toll, but a rough estimate would suggest about a thousand operations [70] . 

There have been several attempts [7T] IT*"] to model the effects of decoherence in 
small NMR quantum information processors. We also used a simple model [""U| to 
describe the decrease in signal intensity from our two qubit implementation of the 
Deutsch algorithm, attributable to relaxation. The experiment is explained in Chapter 
Eland the model itself is outlined in Appendix ??. Interestingly, it is also possible to 
reverse the effects of unwanted system-environment interactions [73*1 IT""] by perturbing 
the system, a technique called dynamical decoupling; or avoiding these interactions 
altogether, by working in decoherence-free subspaces [73]. These noiseless sub-spaces 
have now been experimentally realized in NMR quantum computers [7"""| [77] . 

1.1.11 Gradients 

However, there are sometimes good reasons of intentionally introducing decoherence 
into the ensemble: as the name suggests, the gradient magnetic field [HI [THJ introduces 
a spatially linear magnetic field in the sample. If the gradient strength is G, a point in 
the sample at the spatial coordinate z, will experience a field of Bq + Gz, instead of the 
otherwise homogenous field, B Q , resulting in a position-dependent Larmor frequency, 



If we consider only the spatially dependent part, the two-qubit Hamiltonian during the 
gradient pulse, H g can be written as, 



uj= -j(B + Gz)(l-a). 



(1.33) 




7/Gz(l - oi)l z - j s Gz(l - a s )S z ; 



(1.34) 



and furthermore, assuming a homonuclear system with 07 rs as = cr pa 0, 




jGz(I z + S t ), 



(1.35) 



and the propagator is given by, 



U g (z) = exp (-iHg(z) t) = exp (-i</){z) (I z + S z )), 



(1.36) 
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where <f>(z) = —jGzt, is a space dependent phase, different points in the sample incur- 
ring different phases. Now consider a coherence term \r)(s\ (an element in the density 
matrix), with coherence order p = M r — M s , where M r and M s are the magnetic 
quantum numbers 7 of the states \r) and \s). The propagator transforms this coherence 
to, 

U g \r)(s\Ul 

= exp (-i(j){z){I z + S z )) \r)(s\ exp (i(f)(z)(I z + S z )) 
= exp (— icj)(z)M r ) exp (icj)(z)M s ) \r)(s\ 

— exp (— i(j)(z)p) \t){s\, (1-37) 

indicating that the z-dependent phase also depends on the coherence order: popula- 
tions and homonuclear zero quantum coherences, having p = 0, acquire zero phase, 
irrespective of the spatial position. In the post-gradient measurement, the detection 
coil integrates the coherences originating from across the length of the sample, and if G 
is sufficiently strong, (or the duration of the gradient is long enough,) the exponential 
terms in (jl.37|) will add up to zero (if p ^ 0). The action of an ideal, strong gradient 
field (also called a crush), is therefore to "crush" all the single and double quantum 
terms in the homonuclear density matrix, 

f u □ □ n\ (u ^ 

□ □ □ □ crush □ □ 

^ . (1.38) 

□ □ □ □ □ □ 

□ □ Dj \ Dj 

Gradients used in this way, are often called purge gradients 0. I will sometimes call 
them crush sequences as well. Note that in heteronuclear systems, even the zero quan- 
tum terms will be suppressed. However, researchers have also come up with techniques 
for suppressing homonuclear zero quantum terms; for an example see jTH]- We shall 
use gradient fields in different contexts through the course of this work. 



1.1.12 Quantum state tomography 

Quantum state tomography |H] is the process of determining a density operator p. A 

two qubit p has sixteen elements, and is characterized by fifteen independent param- 
7 For example, the quantum number M r is defined by the operator equation I z \r) — M r \r). 
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eters [80J, for example, the fifteen coefficients in (|1.12|) . and state tomography entails 
determining these parameters. Direct detection of p, gives access to only the eight 
direct observables, (jl.22j) and so, additional experiments are needed, converting the 
unobservable terms into observables. For example, a two qubit p can be characterized 
in four experiments; possible sets are the experiments involving, 



where 1 represents an identity operation, corresponding to direct detection. (The four- 
step tomography also suggests that the nine-step scheme used in [221 12H] is excessive.) 
In some cases, when there is some prior knowledge about the state, fewer tomography 
experiments may also be sufficient. One obvious method is using a crush gradient, 
which forces all single and double quantum coherences to the value zero. In Section 
13.11 I present an example of a one-shot tomography on a state which is constrained to 
have only two independent parameters. Recent demonstrations of state tomography 
in NMR can be found in [5B1 IH1 E2] The reconstruction of density matrices from 
read-outs in an NMR experiment is exemplified in [53] . 



As I have mentioned, NMR is very good at unitary operations, with fidelities approach- 
ing unity. However, its main weakness lies in non-unitary operations, namely initial- 
ization and single-spin read-out. This section discusses the problem of initialization 
and the following section touches upon the issue of measurement. 

1.2.1 Initial states in "thermal" NMR 

Consider a weakly coupled spin system. The equilibrium density matrix of an ensemble 
of n qubits is given by [31 EH HZ] , 



{1,90^,90^,90 I y }, and 
{90I x ,90I y ,90S x ,90S y }, etc., 



1.2 Initialization in NMR 





Peg 



Z 



where, 
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is the partition function and serves for normalizing the state and H is the system 
Hamiltonian. The operator p eq will be diagonal in the basis of the Hamiltonian H |31j . 
(the Zeeman basis). If q, r and s are indices for the energy levels and i is an index for 
the qubits, the individual matrix elements in p eq will be given by, 

(r\p eq \s) = (r\ex V (-H/kT)\s)/j:l n =l (q\ex P (-H/kT)\q) 

6 rs e*jp(-E r /kT) I J2 2 q n =1 exp(-E q /kT), (1.41) 

where E r is the energy of the rth eigenlevel (Zeeman level), 

E r = fuo r and (1.42) 

= YTi=\ m i^i (1-43) 

homonuclear , . 

« u i YZ =x mi= UiM r , (1.44) 

uji = —^iB being the Larmor frequency of the zth qubit, and the spin magnetic quantum 
number jHlj, m« being 1/2 or —1/2 depending 8 on whether the corresponding qubit is 
| a) or In writing ()1.43|) . I assume that the coupling frequencies J (hundreds of 
Hz) are very small as compared to the Larmor frequencies (hundreds of MHz), and in 
writing (jl.44|) . I also neglect effects of chemical shift differences between the qubits. 

In ([1.41)1 . (r\p eg \r) represents the fractional population in the rth level, suggesting 
a Boltzmann distribution of qubits in the 2 n available levels. For a homonuclear sys- 
tem, these levels are distributed in n + 1 "steps" , in a Pascal's triangle fashion. The 
distribution of levels and steps in the energy manifold is exemplified in Table H~2l and il- 
lustrated for two and three qubits in Figure fTTTl The summary is that the distributions 
of levels in steps of a "Pascal's ladder" is a direct result of the addition of the magnetic 
quantum numbers mj. The sum of the magnetic quantum numbers, M r = Yli m i> can 
take up a value from the range, {n/2,n/2 — l,n/2 — 2, ... , — (n/2 — 1), — n/2}. The 
number of steps is n + 1, and if we index the steps k — 0, 1, . . . , n, then the number of 
levels in the kth step is n Cfc, and the M r of the individual levels belonging to this step 
is n/2 — k. 

With this background, we can derive a product operator form for an n qubit initial 

state. This will also be a convenient entry into understanding the problem of initializa- 
8 Without loss of generality, I assume that m = 1/2 for the |a) state and —1/2 for the |/3) state. 
This translates into \a) being a low energy state for positive 7 nuclei such as hydrogen. In fact, 
understanding the signs of frequencies in NMR is a convoluted undertaking; see, for example [41) . 
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no. of qubits no. of levels no. of "steps" levels per step 



12 2 1 

2 4 3 1,2,1 

3 8 4 1,3,3,1 

n 2^ Ti + l l,"Ci,"C 2 ,..., w C ra -i, ra C n 

Table 1.2: Arrangement of energy levels in the nuclear spin energy manifold for n qubit 
homonuclear systems, neglecting coupling interactions. 



M= -1, k=2, "C k =l 

M=0, k=l, "C k =2 

M=l, k=0, "C k =l 

(a) n=2 

M= -3/2, k=3, "C k =l 

M=-l/2,k=2,"C k =3 

M=l/2, k=l, "C k =3 

M=3/2, k=0, "C k =l 

(b)n=3 

Figure 1.1: Conceptual representation of energy levels in steps on a Pascal's ladder, 
and correspondence with the total magnetic quantum number M r . Case (a) is for two 
qubits (n = 2) and case (b) represents three qubits (n = 3). The M r , step index k and 
number of levels in each step, n Ck are written alongside the steps. 
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tion in NMR. Consider a single qubit n = 1; there are two energy levels, each level in a 
separate step, \a) being lower in energy than \(3). Suppose the fractional populations 
in these states are n a and np] from (jl.41|) . n a is given by, 

n a = eyjp(-E a /kT)/(y^ j exp(-E q /kT)) 

q=a,/3 

= exp (— fkOa/kT) /(exp (— Juu a /kT) + exp (— hujp/kT)), (from (|1.42|) ) 

= exp(~hu/2kT)/(exp(-hiu/2kT) + exp(fku/2kT)), (from Oil ). (1.45) 

where u is the Larmor frequency of the first and the only qubit . In the high temper- 
ature limit, kT ^> ftw, n a is approximated as, 

n a » (l - huj/2kT)/{\ - hcu/2kT + 1 + hw/2kT) 

= (l-huj/2kT)/2, (1.46) 

and similarly, 

np*i(l + hjjj/2kT)/2. (1.47) 
With these fractional populations, the equilibrium state can be written as, 

Peg = {K^}} 

1(1- hw/2kT 
2 I 1 + Huj/2kT / 

= i(l 2 -BI,)~I„ (1.48) 

where i3 = fko/kT is what I refer to as the Boltzmann factor. Similarly, for two 
homonuclear qubits, the equilibrium state is 

l(l 4 -B{I z + S z )) ~I Z + S Z , (1.49) 

and for n qubits, 

n n 
i=l i=l 

For a proton frequency of 400 MHz, the high temperature approximation will be 
valid for T ^> 0.02 K, indicating that the equilibrium states (|1.5U|) are an extremely 
good description at temperatures of even a few hundred millikelvins, and therefore, 
must also accurately describe our initial states at room temperature. I shall refer to 
these high temperature equilibrium states as thermal states in the remainder of this 
thesis. 
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Now some important comments on the energy gaps and population differences be- 
tween levels, which are significant because the NMR signal arises from differences in 
level populations. The ball-park figures derived in the following discussion all assume 
a proton frequency of 400 MHz and a temperature of 295 K. The Zeeman splitting 
defined as the energy gap between levels in adjacent steps, is always Hu, resulting in 
the steps being evenly spaced. This splitting is extremely small, a tiny 2.6 x 10~ 25 J, 
and the available thermal energy kT ^4x 10~ 21 J is four orders of magnitude larger, 
resulting in almost equal populations between the levels. Another way of stating this is 
that B is exceedingly small: about 6.5 x 10~ 5 . However, the populations in the adjacent 
levels are not exactly identical, otherwise thermal NMR would not have been possible 
in the fist place. In fact, the fractional population difference between adjacent levels is 
proportional to B. For example, in a single qubit, the difference is B/2 « 3.2 x 10~ 5 : 
which means that for roughly every 3 x 10 4 spins, there will be one extra spin in the 
lower | a) level and if the sample contains 10 15 spins, the net excess of spins will be of 
the order of 10 10 ; not a very large number, keeping in view the very small magnetic mo- 
ment of the nucleus. For all other nuclei (except radioactive tritium, 3 Hi), the Zeeman 
splitting is even smaller. 

1.2.2 Pseudopure initial states 

The foregoing discussion explains the origins of the notoriously low sensitivity of NMR 
arising from the tiny Zeeman splittings: small gaps lead to small population differences 
which in turn lead to small signals. The insensitivity is a problem in analytical NMR, 
especially if the concentration of the sample is small. A lot of effort, therefore, is 
devoted to attaining increased polarization or hyperpolarization and some techniques 
for achieving this will be touched upon in Section 11.51 However, in the context of 
quantum information processing, the low magnetic to thermal energy ratio is also a 
limiting factor in a slightly different way and I now outline this limitation. 

Every computational process must start off by placing the computational register in 
a well-defined initial state; for a quantum computer, we normally expect the (register 
of) qubits to be placed in the energetically ground state |0) = 1 00 . . . 0) = \aa ... a), a 
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pure quantum state. For example, for two qubits, a suitable initial state would be, 

{{1,0,0,0}}= ^ + h+s*+2hs z 

~ I Z + S Z + 2I Z S Z , (1.51) 

having only one non-zero eigenvalue. What nature furnishes instead is the thermal 
state, (jl.49J) ~ I z + S z = {{1,0,0, —1}}, having two non-zero eigenvalues. We imme- 
diately encounter two problems: first, the thermal and desired initial states are not 
identical and second, the amount of the (deviation) thermal state is proportional to 
a very small number, B. How do we then prepare a pure initial state, suitable for 
quantum computing? The answer was provided by Cory et al, [HUES], who suggested 
that it is not necessary to prepare pure states, rather it is sufficient to prepare only 
pseudopure states, that are given by the general form, 

Pps= (l-e)^ + e\^)Mp, (1-52) 

the subscript "ps" and "p" indicating pseudopure and pure states. For a two qubit 
system, a conventional pseudopure state is, 

P PS = (l-e)j+e|00)(00|, (1.53) 

comprising e parts of the desired pure state mixed with 1 — e parts of the maximally 
mixed state. This factor, e, is what I term the polarization, signifying the amount of 
"polarized" pure component in the overall state, the maximally mixed (unpolarized) 
state being NMR silent. The dynamics of p ps are exactly identical to 1 00) (00 1 , and so 
the pseudopure state effectively behaves as a pure state, albeit with smaller polarization 
and hence lower sensitivity. 

How do we prepare these effective pure states? The eigenvalue spectra for the 
thermal and pure states are different, and so some kind of non- unitary operation is 
required [3TJIEE] : as a result, the different approaches for assembling pseudopure states 
all involve some kind of non-unitary averaging of eigenvalues. At the end of the aver- 
aging process, we are left with one non- degenerate eigenvalue. (A pure state has one 
non- zero eigenvalue.) 

The different strategies for pseudopure state generation are reviewed in [T5^ ll9j. In 
summary, the earliest proposals [T71 |H3] , called spatial averaging, employed field gra- 
dients, whereby the averaging is achieved by entangling the spin degrees of freedom 
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with spatial dimensions and later taking a partial trace over space. Another approach 
|H71 15%] . called temporal averaging, performs the averaging by carrying out multiple 
experiments, each experiment using a different initial state unitarily derived from the 
thermal state. This second approach, in a sense, entangles the spin degrees of freedom 
with time and throws away the time information at the end. A third scheme !9Uj. 
called logical labelling, selects a subset of qubits embedded in a higher dimensional 
quantum space, conditioned upon the state of the unselected qubits. This last ap- 
proach is elegant in principle, but the complexity of the pulse sequences has limited its 
applicability to small homonuclear systems. There is yet another efficient strategy j^I] 
for "extracting" an n — 1 qubit pseudopure state from n-multiple quantum coherence 
terms (related to "cat" states [Hj) which are selected using filtration sequences based 
on the use of gradients or phase cycling |3J. 

1.2.3 Climbing "Mount Scalable"? 

A large part of this thesis concerns the preparation of a two qubit NMR based quantum 
computer, so what method do we use to prepare pseudopure states? The answer is, 
simply, we do not prepare pseudopure states, instead our system is directly initialized in 
an almost pure state. A part of the very first DiVincenzo requirement ^T] is scalability 
and the pseudopure strategy does not satisfy this requirement, as I now review. Most of 
the following criticism was anticipated in Warren's paper [§2] and more comprehensive 
reviews on the scalability issue can be found in |701 19*3] , the latter reference being the 
source of the term "Mount Scalable" . 

Warren derived an upper bound on the polarization factor for a pseudopure state 
in a homonuclear system, which I call the Warren bound. The derivation is based on 
the fact that the theoretical maximum polarization achievable from a thermal state is 
upper bounded by the maximum fractional population difference in an n qubit system: 
the maximum difference is between the \aa ... a) and \(3j3 . . . (3) states, which are the 
furthest apart in energy. These "all a" and "all (3" states have magnetic quantum 

numbers n/2 and —n/2, and using (|1.41|) - (|1.44j) . the corresponding difference 9 can be 
9 This difference enables excitation of the multiple quantum coherence of order n 4 . In reference 
|B2| . the expression (jl.54|) does not carry the negative sign, however, this sign is consistent with our 
assumption that for positive 7 (and negative ui) nuclei, \a) is the lower energy state, resulting in 
n a u a > n a ii f}- 
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derived as, 

^aii a — n n\\ p — ( ex P ( — nhjj/2kT) — exp (nhu / 2kT)*j / Z 

= -2smh(nhu/2kT)/Z. (1.54) 

The partition function Z can be written in compact form as the Binomial expansion, 

n 

Z = ^ n C k (exp(-huj/2kT)) n ~ k (exv(hu/2kT)) k 

k=0 

= (exp(-huj/2kT) + exp(frLu/2kT)) n 

= 2 n cosh n {hw/2kT), (1.55) 
and therefore, Warren's bound is given by, 

Kn a ~ p) jZ = -2 sinh (nfru/2kT) / (2 n cosh" (Hu/2kT)) , (1.56) 
which in the limit of high temperature and ignoring the negative sign, is, 

71 71 

¥ (hu/kT) = —B, (1.57) 

where I have used the limiting values of the hyperbolic functions, linx^osinhx = x 
and lim^o cosh x = 1. 

The exact polarization extracted from the the thermal state, depends on the method 
used to prepare the pseudopure state, however in all cases, the Warren bound (jl.57j) 
sets an upper limit. In the high temperature limit, the bound results in an exponential 
decrease in polarization as we linearly increase the number of qubits; and therefore to 
defeat the exponential loss in sensitivity, exponentially more resources are required, 
such as exponentially larger samples. This is not only impractical but also defeats 
the fundamental gain expected (and desired) from a quantum computer, that is ex- 
ponentially bigger Hilbert dimensionality with at the most, polynomial resources. It 
appears that thermal NMR based on assembling pseudopure states is not scalable and 
this statement is the gist of the initialization problem. 

1.2.4 Concentrating polarization 

However, the literature suggests that the problem of low polarization and scalability 
can be sidestepped, using a computational (technology independent) approach outlined 
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by Schulman and Vazirani |94j . The method centres around polarization concentra- 
tion: given an ensemble of n weakly polarized qubits, it is possible to extract k < n 
hyperpolarized qubits, at the expense of maximally randomizing the remaining n — k 
qubits. The scheme for extracting pure qubits can be represented as, 

®((l-e)y + e|0)<0|)-^ ((g)|0>(0|)®((g)^), (1.58) 

n k n—k 

denoting the tensor product between j qubits. The algorithm operates on classical 
states and conserves the classical entropy S [8 of the system. From this condition, an 
upper limit on the number of pure qubits extracted can be determined: the entropy of 
a single qubit with polarization e (and fractional populations (1 ± s)/2) is 

S(s) = -(^±£) lg 2 (^) - (^) lg 2 (^), (1.59) 
and as the entropy is conserved, we can write, 

nS(e) = kS{l) + {n-k)S(0) (1.60) 

= fc(0) + (n-Jfc)(l) 
=>- k = n(l-S(e)) 

= n(l + (^)lg 2 (^) + (^)lg 2 (^)), (1.61) 

(1.62) 

and using the expansion lg 2 (l ± x)/2 = — 1 ± xj In 2 — x 2 / (2 In 2) ± 0(x 3 ) for small x, 

TIE 2 

k « — — - « ne 2 (for large n). (1.63) 

The equation indicates that we can extract k pure qubits from k/e 2 low-bias qubits; 
the scaling is only linear in k. This may seem to be an exciting possibility for NMR 
but taking a more practical point of view, the restrictively small values of e ~ 10~ 5 in 
thermal NMR, require a prohibitively large (~ 10 10 ) number of initial qubits, even to 
extract a single pure qubit. However, despite this limitation, an NMR demonstration 
of the basic building block in the polarization boosting has already been demonstrated 

m 

Another variant of the above technique is the so-called algorithmic cooling method 
|9"6~] which is an entropy non-conserving algorithm for polarization concentration. The 
method beats the Shannon bound ()1.6()|) through the use of rapidly relaxing qubits in an 
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open quantum system and an initial implementation of this non-adiabatic version has 
also been reported [HZj- However, it seems that both of these computational schemes 
will become practically useful in small sized spin systems, only when used in conjunction 
with higher initial polarizations. 

1.2.5 Re- initialization 

Quantum error correction is an attempt to counteract the effects of decoherence, 

reversing errors once they have occured. Thermal NMR is not well suited for this pur- 
pose, because of two major reasons: first, most error detection and correction schemes 
work by encoding a physical qubit in several logical qubits, the increased number of 
qubits reducing the sensitivity; and second, repeated error correction requires one to 
reuse the ancillary qubits and currently, no practical schemes are known for resetting 
NMR qubits into pure states at will. Thermal NMR cannot realize pure states initially, 
let alone provide a continuous supply of pure qubits during the experiment and so, even 
though, simple error detection [100J and correction |lUlj have now been implemented 
in small NMR systems, they mostly remain proof of principle experiments. 

1.3 Measurement in NMR 

Liquid state NMR does not implement projective measurements jHj on single qubits. A 
projective ("strong") measurement on a qubit in the state (co|0) + Ci|1))(cq(0| + (1 1 ) , 
yields |0) with probability |co| 2 and |1) with probability |ci| 2 , leaving the state in |0)(0| 
or |1)(1|, also showing that such measurements constitute one (though not the only) 
way of initializing the qubit. In liquid state NMR, we do not have access to single qubits 
but rather ensembles of qubits, and the so-called "weak" measurements performed in 
these systems provide ensemble averages \TJ\: the final result is a noisy expectation 
value of a spin operator, such as (I x ) or (I y ), taken over all molecules in the sensitive 
region of the detection coil. Furthermore, the wavefunction does not collapse into one 
post-measurement state as the system is only feebly affected by the detection process 
|1T]I1(J2| . In short, it appears that these weak measurements are statistically equivalent 
to standard, single-qubit measurements, and will suffice, until we make enough progress 
to start actively considering resetting qubits for quantum error correction. 
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It is not out of place to comment upon a few research initiatives aimed at single-spin 
or single-molecule detection. Single sub-atomic particles have small magnetic moments 
pH] and are difficult to measure in normal circumstances. There has, however, been 
recent progress in detection of single electron spins using the spin-dependent magnetic 
force on a ferromagnetic tip |103j and spin-dependent charge transport |lU4j . Direct 
spin detection of protons and nuclear spins is even more difficult, as their gyromagnetic 
ratios are much smaller. The detection of single nuclear spins, therefore, requires 
some kind of "amplification" of nuclear Zeeman splittings by correlating them with 
the wider energy gaps present in another degree of freedom. For example, in the first 
optically-detected magnetic resonance experiments |105| HU6j , the magnetic transitions 
(microwave) in the electronic spectrum of a single molecule were indirectly detected as 
optical transitions, which are intense and can be monitored with much higher quantum 
efficiencies. Another interesting experiment in this regard is the optical detection of 
nuclear magnetic transitions in a single quantum well |lU7j . It is widely felt that any 
envisaged single-nuclear spin detection scheme is likely to involve detection of optical 
photons. 

I have now discussed liquid state NMR in terms of the DiVincenzo criteria. There 
remains another important and more fundamental aspect of thermal NMR that has 
even raised questions about its acceptance as a truly "quantum" process. In the next 
section, I take up the widely debated issue of entanglement in liquid state NMR. 

I. 4 Entanglement in NMR 

Entanglement jS] is undoubtedly an important resource in quantum information pro- 
cessing; some researchers would argue that it is the "characteristic trait of quantum 
mechanics" 108 . Entanglement is interesting because it enables quantum states to 
exhibit correlations that have no classical analog. Of its many example applications in 
quantum information science, entanglement enables quantum teleportation 1 1 001 1110] . 
superdense coding |111) 1112) and secure communication |1131 1114j . 

A bipartite pure state will be entangled if its Schmidt number is greater than one 
jS] and a mixed state of two or more qubits will be entangled if it cannot be written 
as the convex sum of the direct products of the component density operators (Section 

II. 4.1)1 . Several operational criteria for detecting and quantifying entanglement will be 
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taken up in the following subsections. 

As an example of the consequences of entanglement, consider the maximally entan- 
gled state, 

l0 + )= ^(|01> + |10». (1.64) 

Tracing out the second qubit, the reduced density matrix jH] of the first qubit will 
be the maximally mixed state, pa = 1/2 {{1, 1}}, implying that measuring this qubit 
along any axis would result in a random result, and thus measuring one subsystem 
alone provides no information about the preparation of the composite state, (jl.64|) . In 
this sense, the entangled state contains "hidden information" which cannot be revealed 
by measurements on either qubit. However, if both subsystems are measured and the 
results are compared, they will be perfectly correlated. 

Conventional NMR based QIP uses highly mixed, separable (non-entangled) states. 
In this section we explore the relations between entanglement and the purported power 
of quantum computers, more specifically, NMR quantum computers. For example, can 
NMR devices exhibit entanglement? And are current devices truly quantum? In fact, 
most of the work in the current thesis addresses the two-fold problem of initializa- 
tion and entanglement. Subsequent chapters will discuss in detail, our experimental 
implementations, of initializing a two qubit system in an entangled state. A gentle 
introduction to the role of entanglement in NMR can also be found in |115j . 

1.4.1 Separable and entangled states 

Traditional QIP experiments using liquid state NMR furnish and employ separable 
states. A two qubit pure state will be separable if it can be written as a tensor product 
[HI EI] °f the component pure states, 

\x)ab= |^}a®|0)b, (1-65) 

where A and B label the qubits. Examples of separable states include 1 00) = |0) <g> |0) 
and (|01) + \ll))/y/2 = ((|0) + \l))/y/2) <g> |1). A pure state is non-separable, or 
entangled, if it cannot be brought into the form ([1.650 ; examples are the Bell states 
[H], defined as, 

|^)= (|00}±|11))/V2 (1.66) 
|^)= (|01)±|10»/V2. (1.67) 
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However, as pointed out earlier, states encountered in traditional NMR are mixed 
and have exceedingly small polarizations. Separability for such mixed states is, at least, 
well- defined: a two qubit mixed state p will be separable if it can be decomposed as a 
weighted sum of product states, 

P = ^Pi{pA® Pb) v with pi > and ^Pi = 1. (1.68) 

i i 

The product state (p^ ® Pb) { , prepended with a non- negative classical probability Pi, 
has the property that the A and B quantum sub-systems are independent: knowledge 
of one sub-system does not depend on or affect the other sub-system. Product states 
are not entangled; and neither are their convex combinations, such as (J1.68)) . All states 
in low-polarization NMR can be written in the form (J1.68)) . and are therefore, always 
separable. 

A state we shall frequently encounter in this work is the Werner state |116j . For 
two qubits the state is given by, 

Pe = (1 - + £\^ent) v0ent|, (1.69) 

where \ip e nt) is one of the four maximally entangled Bell states, (jl.66|) and ([1.67)1 . 
Whether p e is entangled or not, depends on the precise value of e, and is investigated 
in Section 11.4.41 The headline result is that the Werner state (11.69)1 will be entangled 
when e > 1/3. 

1.4.2 Separability bounds for mixed states in NMR 

Until now, I have defined separability and given examples of two qubit entangled states. 

Can we also come up with general separability tests for n qubit states, where n > 2? 

This question is important in understanding the scalability characteristics of quantum 

computers, and in context of NMR, was first addressed by Braunstein et al. [1 17j . The 

major results of their work and some interesting observations are summarized in this 

section 10 . A general method 11 for deriving the lower bound and constructing explicit 

decompositions of separable states can also be found in [120J. 

10 It should be noted that since Braunstein's original work, tighter bounds have been proposed 

EBI Eg. These new bounds modify the finer details, but do not affect the overall broad results. 
11 The method is based on decomposing a density operator in an over-complete basis, and based 

on the constraints on the operator eigenvalues, looking for conditions that guarantee non- negative 

probabilities in (|1.68J) . 
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NMR deals with highly mixed states. The maximally mixed state M n = 1/2" is 
clearly separable, as are states lying sufficiently close to M n . Braunstein and co-workers 
determined |117j lower and upper bounds on the size of the separable neighbourhoods 
of M n , for arbitrary n. The concept is illustrated in Figure 11.21 indicating that the 
bounds result in three distinct regions: S, where all states are provably separable, 
E, where entanglement can definitely exist, or in other words, where the states are 
unitarily equivalent |121j to entangled states, and a region in between, ES, where it 
is not clear whether entangled states can or cannot exist. Consider the n qubit state, 
with a "distance" e from M n , 

p= (l-e)M n + e Pl , (1.70) 
pi being an arbitrary density matrix. For sufficiently small e, 



I large n 2 

1 _|_ ~ 4 T 



e<ei = — « — , 1.71 

— ' 1 i o2n-l /In' v > 



p lies in the region S. Similarly, for sufficiently large values 12 of e, 

1 large n 1 

£ > £ u= j - 2 „ /2 ~ 2^2' I 1 ' 72 ) 

entangled states can exist. These bounds are shown in Figure ESI 

Some important observations emerging from the Braunstein and Warren bounds 
are the following. 

1. As the number of qubits increases, the size of S shrinks, increasing the likelihood 
of entanglement, even for reasonably small values of e. 

2. For a typical value of B ~ 10~ 5 , and the modest number of qubits presently 
employed in NMR quantum computers, (the current record is at seven qubits 
|2*U].) all states lie in the region S. This has lead to interesting discussions raising 
suspicions about the very nature of liquid state NMR as a truly "quantum", 
computing paradigm. 

3. Fixing £>, we can estimate the number of qubits at which the S to ES cross-over 
will take place. For B = 10~ 5 , this transition would occur with fourteen qubits. 



12 The upper bound is determined by projecting an arbitrary 2 d dimensional state of two qudits 
onto a four dimensional Werner state of two qubits, and using Werner's separability criterion |122j to 
estimate the threshold. 
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Figure 1.2: Lower e\ and upper e u bounds as constructed by Braunstein and co-workers 
|117j . The regions S, E and ES indicate the existence of definitely separable, definitely 
entangled and "maybe-entangled" states. The maximally mixed M n lives in the middle 
of S. The continuous S-ES boundary belongs to S, and the dashed ES-E boundary 
belongs to ES. 

4. If instead of fixing B, we choose to fix n at seven qubits, then B should be at 
least 0.002 for the transition into ES. On modern high-field spectrometers, this 
translates into a maximum temperature of about 10 K; most liquid samples would 
be frozen at this temperature. 

5. Figure IT731 also suggests that with low polarization NMR, it will never be possible 
to cross over into the provably entangled region E: a radically different approach 
is therefore needed to realize entanglement in liquid state NMR. 

1.4.3 (Quantum) nature of liquid state NMR (quantum) in- 
formation processing 

The above findings have led to considerable concern about the nature of liquid state 
NMR: the statistics of experimental outcomes from separable states can be described 
perfectly well using classical probability distributions, and therefore, in its very essence, 
liquid state NMR may not be quantum at all |123j ! This might seem to be the end of 
the story for liquid state NMR quantum computing, but luckily, the details are much 
subtler. 
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Figure 1.3: Lower (ej) and upper (e u ) separability and Warren bounds as a function 
of the number of qubits n. The ordinate scale is logarithmic. The line with the open 
circles indicates £/, the S-ES boundary and the one with filled circles represents e u , 
separating ES and E. The line identified with squares is the Warren bound, plotted 
assuming B = 1CT 5 . 
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To start with, there are differing "definitional" stances on the "quantumness" of 
liquid state NMR, views revolving around the question of whether entangled states are 
a prima facie signature of quantum behaviour, or if the real quantum character lies in 
quantum mechanical dynamics, or possibly somewhere else. No doubt, the quantum 
states are separable, but what about the dynamics? Can quantum information also 
"reside" in states that are not entangled, as discussed in the context of quantum discord 
[124J? Finally, can we come up with efficient classical models to simulate the dynamics 
of separable states? 

In short, the role of entanglement in quantum computing per se, remains an open 
question |115t 1125] . At the one end, entanglement has been viewed as the source of the 
enhanced information processing capabilities of quantum computers, a view resonating 
in |126j . But at the same time, liquid state NMR (quantum) information processors 
routinely implement quantum dynamics governed by the Liouville-von Neumann Equa- 
tion |2], even though they do not exhibit entangled states. Having said this, we must 
also remember that efficiency arguments must be analyzed in light of the fact that the 
enhanced power of quantum computers is only an assumption, although a very tempt- 
ing one, and except for the black box or oracle case, there is no proof as yet of the 
increased power of quantum processors. Second, quantum algorithms outperform only 
the best known classical algorithms. For example, in the case of the factoring problem, 
one day somebody might also come up with a classical algorithm for efficient factor- 
ing |127j . In this case, most of our fundamental perspectives on quantum information 
processing would be subject to major change anyway. In light of this discussion, we 
can safely conclude, that the essential nature of high temperature liquid state NMR is 
an open-ended question, and is intertwined with the supposed role of entanglement in 
conferring arguably better information processing capabilities to quantum computers. 

The present work, however, by-passes this interesting discussion by preparing de 
facto entangled states; the first demonstration of entanglement in liquid state NMR. 
As a result, most of the charges made against the very nature of liquid state NMR 
implementations can no longer be directed against our two qubit quantum computer. 
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1.4.4 Detecting entanglement in two qubit states 

One of the main claims made in the forthcoming chapters is that we achieve entangled 
states in liquid state NMR. In this respect, one might ask the following questions. How 
can entanglement be detected in our two qubit system? Finally, can we go one step 
beyond detection, and actually quantify the amount of entanglement? This section 
deals with these questions. Useful primers on the subject of detecting and quantifying 
entanglement are the references |128t [129]. 

Positivity of partial transpose 

For a two qubit state, entanglement can be detected, for example, through the so- 
called positivity of partial transpose (PPT) test, first proposed and exemplified in 
|122j . Consider the separable mixed state (jl.68[) . Taking the transpose of the reduced 
density matrix pH] . pa, results in an overall state, 

P Ta = (pa ® Pb) v with pi > and >^Pi = 1, (1.73) 

i i 

where T denotes the transpose operation. As the reduced matrix pa is hermitian (being 
a valid density operator [H]), p T A = p* A is also an hermitian, unit trace matrix and hence 
another legitimate density operator. This, in turn, implies that the state p TA is also 
a valid density operator and has non-negative eigenvalues. The PPT test, therefore, 
detects non-separability by searching for negative numbers in the eigenvalue spectrum 
of the partial transpose 13 of p. 

The PPT is a simple operational test and I illustrate it with a relevant example. 
Consider the Werner state, 

p £ = {l-e)^+e\i,-)^-\ 
( ±f ^ 

u 4 2 u 

o -§ if o 
v o 

13 It is equally valid to find p TB instead of p TA . 
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(1.74) 
(1.75) 



and taking the partial transpose with respect to the first qubit, 
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we obtain an eigenvalue spectrum, {(1 + e)/4, (1 + e)/4, (1 + e)/4, (1 — 3e)/4}. For le- 
gitimate values of e, the first three eigenvalues are always positive, whereas the fourth 
eigenvalue will be non- negative only He < 1/3. This immediately suggests that an 
inseparable or entangled p £ will always have e > 1/3. In fact, the 1/3 threshold can 
also be deduced by setting n = 2 in Braunstein's upper bound, (jl.72|) . Horodecki, 
Horodecki and Horodecki went on to show [130J that the positivity of the partial trans- 
pose is not only a necessary, but also a sufficient requirement for separability of systems 
of dimensions 2x2 and 2x3, but the result cannot be extended to higher dimensional 
systems. (For example, there exist |131j entangled states which satisfy the PPT test; 
such states are called bound entangled states, and are undistillable |132j .) Multi-partite 
entanglement, is in fact, an open theoretical challenge in present-day quantum infor- 
mation science, and many general results are still unknown or are being investigated 
(for example, see |133j ). However, for two qubits (a 2 x 2 system), the situation is 
rosier, and luckily, we shall be primarily concerned with two qubit (four dimensional) 
quantum systems in most parts of this thesis. 

1.4.5 Quantifying entanglement in two qubit systems 

As indicated earlier, qualifying entanglement is an arduous task, especially for higher 
dimensional mixed systems. How can we expect the quantification to be any easier? 
Several conditions for a "good" entanglement measure have been compiled |8Ul29 |IT3*4*j . 
and there do exist a range of measures, satisfying some or all of these conditions. Once 
again, for two qubit systems, the situation is relatively simple and well understood. 

Concurrence and entanglement of formation 

According to the PPT criterion, negative eigenvalues of p TA are a signature of en- 
tanglement; the negativity of the most negative eigenvalue is one straightforward mea- 
sure. However, for another entanglement measure, called the entanglement of formation 

34 



(EOF) |135j . there also exists a closed form mathematical expression pf3o] (in the case 
of two qubits). We denote this quantity as Ef, and is our measure of choice for detect- 
ing and quantifying entanglement: a non-zero Ef indicating the existence of entangled 
states. A prescription for finding Ef |136| for a state p is outlined below. 

1. Determine the "spin-flipped" operator, p = (cr y ® a y )p*(a y <g> a y ), where p* is the 
complex conjugate of p. 

2. Find the non-Hermitian operator, p' = pp. 

3. Compute the positive square roots of the eigenvalues of p' and sort them in 
decreasing order, {Ai, A2, A3, A4}. 

4. Find the concurrence defined as C = maximum(0, Ai — A2 — A3 — A4). 

5. The value of C is then used to calculate the value of 

x= 1(1 + 713^2). ( L77 ) 

6. We then use x to calculate the entanglement of formation, through the formula, 

E = -x lg 2 x - (1 - x) lg 2 (1 - x), (1.78) 

reminding us of the classical binary entropy [E] function. Note that in computing 
(jl.78)) . we use the convention lg 2 = 0. 

The concurrence and entanglement of formation Ef, can both be used as legitimate 
measures of entanglement. Both measures are zero for a separable state and non-zero 
for entangled states. For maximally entangled states, Ef is 1 and C is 1. Furthermore, 
Ef is a monotonically increasing convex function of C [136], the functional dependence 
is shown in Figure H~H 

1.4.6 Entanglement in singlet-triplet mixtures 

In this section, I determine Ef for a special class of two qubit density matrices: mixtures 
involving four states |^")(^~| = So, |^ + )(^ + | = T , |00)(00| = T x and |11)(11| = T_i. 
The first two (So and To) are the Bell states 1)1.67)1 and the other two (T\ and TLi) are 
separable states. So is conventionally called the singlet and T , T±i, the triplets. I also 
define T m as, 

T m = + T_i), (1.79) 
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Figure 1.4: The entanglement of formation Ef as a function of concurrence C. 

and observe that T m is also an equal mixture of the Bell states and 

These S/T mixtures have a special place in the present work, as our initialization 
experiment, to be discussed in Chapter El prepares states containing combinations of 
only 5*o, T and T m . The ideal experiment prepares So, but the triplets mix in, in 
varying proportions, resulting in an S/T mixture. Whether the resulting mixture is 
entangled or not, depends on what is mixed in and in what amounts. Based on Ef 
measurements, we can outline some important results. These results are helpful in 
determining the correct threshold for entanglement. 

1. The Ej for the pure So is obviously 1. 

2. Now consider the "binary" So /To mixture, eSo + (1 — e)T , where s < 1. The 
resulting Ef as a function of the amount of singlet, e, is shown in Figure IT3T a). 
showing that Ef is always non-zero, except when the the amounts of So and To are 
perfectly balanced, i.e., when e — 1/2. Interestingly, adding an entangled state 
(T ) to another entangled state (So) reduces the overall amount of entanglement. 
This can also be seen in light of the convexity property |129j of the entanglement 
of formation, 



E f (eS + (1 - e)T ) < eE f (S ) + (1 - e)E f (T ), for e < 1. 



(1.80) 
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The take- home lesson is that the Sq/T mixture will always be entangled if e > 
1/2- 

3. Figure IToT b) shows the Ef for an Sq/T\ mixture: like the case for T , adding T\ 
also reduces the entanglement but it is not as effective in "pushing" the state to 
the separable region, and the mixture is entangled for e > 0. 

4. If we mix in T m , the state is entangled for e > 1/2; the Ef is plotted in Figure 
11.5( c). This indicates that T m , a balanced mixture of separable states, is more 
effective in destroying the entanglement of the singlet, as compared to adding T\ 
or T_i on its own. In fact, T and T m are the most effective states in destroying 
the entanglement in the singlet state. The motivation for this argument is also 
discussed in Appendix ??, Section ??. 

5. Finally, the "ternary" S /T /T m mixture, 

p = eS + (1 - e) ( X T + (1 - x)T m ) , where < e, X < 1, (1.81) 

will always be entangled for e > 1/2, irrespective of the value of x- This is 
illustrated in the two-dimensional plot, Figure PToT d). which also points towards 
additional subtleties, as now, e > 1/2 is not a necessary but only a sufficient 
condition. This is suggested by the non-zero Ef for large values of Xi even when 
e < 1/2, indicating that this time, the mixture is "reflecting" the entanglement 
present in the T state. 

6. For arbitrary mixtures, the state is always entangled when e > 1/2; Ef depends 
only on e, and is independent of the relative proportions of T and T m . A simple 
proof to this hypothesis is presented in Appendix ??, Section ??. 

7. I now come back to the Werner state, ([1.74)1 . for which I have already shown that 
p e will be entangled for e > 1/3. In fact, the Werner state can also be viewed as 
a ternary mixture, by decomposing the maximally mixed state into the singlet 
and triplets, 

I4 Sq + T + 2T m 

T " 1 ' (L82) 

enabling us to re-write p £ as, 

p £ = (If )(S + T + 2T m )+eS 

= (if )(T + 2T m ) + (^)S , (1.83) 
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Figure 1.5: Entanglement of formation Ef for singlet triplet mixtures: (a) So and To, 
(b) 5*o and T 1; (c) So and T m and (d) S , T and T m . The Ef versus e plot for an 
S /T_i mixture is identical to (b). 

where the co-efficient of Sq, signifying the amount of singlet is now (3e + l)/4, 
instead of e. The e > 1/3 condition deduced above, then translates to the 
condition (3e + l)/4 > 1/2, which is identical to the threshold derived above 
for S /T /T m mixtures. Thus, states involving 5*0, T and T m , (including the 
Werner states), will always be entangled if the amount of the singlet exceeds the 
entanglement threshold of 1/2. (A formal term for the amount of singlet in a 
state p is the singlet fraction |135j and is defined as, 

F= (1-84) 

A more concise proof of these arguments can also be found in |137j : the reduced 
density matrices [31] for the state p, (jl.81|) . are the maximally mixed states, I2. It has 
been shown that a state p with maximally disordered subsystems, will be separable if 
and only if its largest eigenvalue is at the most, 1/2. The state p has an eigenvalue /, 
ensuring inseparability for f > 1/2. 

1.5 Methods for enhancing polarization 

It is clear that pure initial states, also lying above the entanglement threshold cannot 
be achieved without resorting to some means for enhancing the polarization. Several 
methods have been proposed and implemented for reaching non-Boltzmann spin popu- 
lations, and most of them have been motivated by the need for enhanced signal-to-noise 
ratios in conventional NMR. I now briefly review some of these methods in the context 
of quantum information processing. 
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1.5.1 High fields and low temperatures 

One of the most obvious methods for increasing the polarization, achieving almost 
pure states, would be to increase B = tko/kT. Calculations with the single qubit 
fractional populations at room temperature, suggest that to achieve a purity of 99% 
(n a = 0.99 and rip = 0.01), a magnetic field of roughly 200000 T is needed which is 
clearly impossible to achieve in the foreseeable future. Alternatively, using a reasonable 
field of 9.5 T requires a temperature as small as 10 mT, which is technologically possible, 
but no liquid exists at this temperature! Therefore, it seems that the naive approach 
is not feasible and a more indirect method is needed. 

1.5.2 Polarization transfer from nuclear spins 

One method to increase the polarization of a qubit would be to "borrow" it from 
another more polarized qubit. Polarization transfer techniques exploit the small pop- 
ulation differences in thermal states and the gains in sensitivity are proportional to 
li/lS) the gyromagnetic ratios of the sensitive and insensitive nuclei. These tech- 
niques, for example, involve the incoherent nuclear Overhauser effect HI U\ or the 
coherent INEPT (insensitive nuclei enhanced by polarization transfer) sequence and 
its refocussed version, jSHHEI; the latter two effectively acting as a kind of SWAP 138 
between the sensitive and insensitive nuclei. However, in these experiments, not only 
are the gains small for most nuclei of interest, but also, it is not clear how they can 
be of any use in enhancing polarizations in X H nuclei, which already have the largest 
gyromagnetic ratios (besides 3 H). Due to these very limitations, these methods are 
likely to be of little value in realizing hyperpolarization for NMR quantum computers. 

1.5.3 Polarization transfer from electron spins 

Instead, it is more promising to make use of the hyperfine interaction: highly polarized 
electrons cannot only facilitate in indirectly detecting nuclear spins (Section II .3|) . but 
can also contribute in the converse manner by transferring polarizations to the nuclei. 
These techniques have been given the generic name of dynamical nuclear polarization 
(DNP) PP, and are closely related to the original Overhauser effect [139, 1140] : one of 
the spins in a coupled pair is perturbed from equilibrium and cross-relaxation effects 
[3J 0] result in enhanced polarization for the unperturbed spin. In DNP, the pair of 
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spins are a nucleus and an unpaired electron. Another related experiment is ENDOR 
(electron nuclear double resonance) [141 j . which traditionally involves transferring 
electron polarizations to the nuclei, performing an NMR experiment and detecting the 
NMR transitions indirectly, by transferring nuclear polarizations back to the electrons. 
For example, using the ENDOR approach, Mehring et al. recently demonstrated jH2] 
preparation of pseudoentangled states between nuclear and electron spins. However, 
all of these methods are mostly employed in the solid state and their real potential in 
the liquid phase is only recently being explored; for example, it has been demonstrated 
that during rapid dissolution of a solid sample polarized by the DNP effect, the nuclear 
polarization is preserved, and can result in enhancements of four orders of magnitude 

ma- 

Optical pumping and spin exchange 

Another important method that produces non-Boltzmann spin distributions is optical 
pumping followed by spin exchange |143| I144j : circularly polarized laser light is used to 
pump electrons in alkali metal vapour (such as Rb) into a non-equilibrium distribution, 
increasing the electronic spin polarization by as much as four orders of magnitude. The 
electronic polarization is then transferred to the nuclear spins in a noble gas such as 
129 Xe via the hyperfine interaction. Optical pumping is, therefore, another example 
of transferring polarization from electrons to nuclei. The method, however, is limited 
in its applicability to quantum computing, because only a few noble gas isotopes are 
known to be polarizable by optical pumping, and being chemically inert, they do not 
normally participate in bond formation, resulting in isolated systems unsuitable for 
universal quantum information processing. In principle though, the spin polarization of 
the Xe nuclei can be transferred to protons using SPINOE (spin-polarization induced 
nuclear polarization) |145t I146| I147j — without the need for irradiating the spins, as 
in a conventional Overhauser experiment — but the transfer efficiency remains very low 
(nearly 10%). One demonstration of using optical pumping to obtain a ten-fold increase 
in polarization over the thermal state, and subsequently implementing a quantum 
algorithm is given in [148] . 
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1.5.4 Chemically induced dynamical nuclear polarization 

CIDNP (chemically induced dynamical nuclear polarization) [3J |3] is a misleading 
name for another technique for generating enhanced magnetic ordering; the method is 
unrelated to the dynamical nuclear polarization discussed above. The effect is based 
on the experimental observation that nuclear spins can "steer" the course of a chemical 
reaction; this seems quite a remarkable fact because the nuclear spin Zeeman energies 
are negligible as compared to bond energies, yet the former can exert a decisive influence 
on bond formation. I present a simple description of this phenomenon. Suppose we are 
given a pair of free radicals A* and B*, where * represents an unpaired electron and 
B is a spin-1/2 nucleus, that can exist in the state |a) or \f3). Also suppose that the 
radicals are "caged" together in solution, and the pair of lone electrons are in a triplet 
state (Tq,T±i). Recombination inside the cage (geminate recombination) takes place 
only if the pair of electrons can be brought into the singlet state (So), as required in a 
stable, bonding molecular orbital |149j . The reaction, 



therefore, can happen only if there is a concomitant triplet to singlet conversion of the 
electronic states. How can such a conversion take place? The nucleus B is hyperfme 
coupled with its electron and causes the electronic spin state to precess at a certain 
frequency, allowing the parallel T configuration to evolve towards the anti-parallel So, 
and the conversion will be more favoured if B is in one particular spin state rather than 
the other. Suppose the \a) state encourages the T — > So conversion. This implies that 
most geminally recombined radicals, that is the molecules AB, will have the B nuclear 
spins predominantly in the \a) state. Likewise, the remaining molecules will eventually 
escape the cage and will recombine instead with the solvent or some other species, 
possibly forming another molecule CB, this time with preferential magnetic ordering 
in the state. At the end, we have two distinct manifolds comprising spin-polarized 
molecules. For QIP, the prospects of CIDNP are as yet unexplored. For example QIP 
requires multiple polarized spins in a single molecule and it is not clear how this can 
be realized by CIDNP. 



A* + B* 



AB 
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1.5.5 Para-hydrogen derived hyperpolarization 

The para-hydrogen induced polarization (PHIP) effect has been used not only to obtain 
large enhancements in spectroscopy |15Uj and imaging |151j applications, but has also 
demonstrated initial success in quantum information processing j^Hl EH 11521 1153| 1154] 
in the liquid state. The following three chapters of this thesis will describe the PHIP 
effect (Chapter EJ); employing it to prepare pure, entangled states (Chapter EJ); and the 
first demonstrations of quantum algorithms on the pure state NMR quantum computer 
(Chapter ij. 
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Chapter 2 

Para-hydrogen induced 
polarization 

". . . NMR is the premier spectroscopic example, not of quantum mechanics, but of 
quantum statistical mechanics • • • •" The fact that most high temperature NMR 
quantum (statistical) computers, are in fact, "ensemble processors" with spins 
statistically distributed in closely spaced energy levels, requires the generation of pseu- 
dopure states. However, we have also seen that the pseudopure approach cannot lead 
to computers of useful sizes. It appears that quantum statistical mechanics limits the 
usefulness of NMR as a viable technology for future computing, and the opening remark 
to this chapter, is a reflection of this view. 

Ironically, quantum statistical mechanics also helps in overcoming some of the ob- 
stacles it creates for itself in the first place: the connection between spin statistics 
and sensitivity enhancement in NMR |155j . was known even before NMR celebrated 
its debut in the quantum computing arena. One possible solution to enhanced NMR 
sensitivity exploits the coupling between nuclear spin and molecular rotation and the 
coupling obeys rules laid out by quantum statistical mechanics. The origin of the 
hyperpolarization lies in the almost perfect spin order, sometimes also called the sym- 
metrization order |156J of the nuclear spin state derived from para-hydrogen. Our 
motivation in using this spin order is in producing nearly pure states in liquid NMR, 
thereby initializing the quantum computer in a well-defined state and realizing entan- 
glement in high temperature liquids. 

This chapter is a foray into the fundamentals of the spin order which arises as 
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a result of the strict correlation between the nuclear spin and molecular rotational 
degrees of freedom. I shall explore the symmetry requirements on the wavefunction 
of a diatomic molecule, such as dihydrogen. These requirements apparently restrict 
"homonuclear freedom" , but in effect, cool the pair of nuclear spins into a pure quantum 
state. The symmetrization postulate j!57t I158| I159j results in the existence of two 
distinct spin isomers of the molecule; this effect is discussed in Section |2~T1 Sections 12.21 
and 12.31 explain the relative proportions of the spin isomers and outline some methods 
of isolating the para isomer. For quantum computing, the exceptionally high purity 
"locked" in the para-hydrogen molecule has to be transferred to a suitable molecule in 
solution; the underlying chemistry is described in Section l2~4l A thorough description 
of the nuclear spin states in para-hydrogen and its derivatives in solution requires a 
formal treatment in the language of density operators. This formalism is spelled out 
in Section The hydrogenation can take place outside or inside the strong magnetic 
field leading to different descriptions. These experiments are outlined in Sections 12.61 
and 12.71 The signal enhancements in these experiments are described in Section 12.81 
Depending on the details of the hydrogenation and detection, variants of the in situ 
experiment are possible, which are explained in Sections l2.10H2~T2l Finally, I present a 
brief summary in Section ^. 131 In short, this chapter prepares the necessary theoretical 
background for the initialization experiment to be discussed in Chapter El 

One of the most comprehensive reviews on para-hydrogen induced polarization 
(which I shall define later) is |!6Uj . while one of the earliest is |161j . Another review 
focusing on the chemical applications is |15Uj . 

2.1 Spin isomers of hydrogen 

In common usage, the word "hydrogen" can refer to either the simplest atom, H, or 
the diatomic molecule, H 2 . To keep the following discussion clear, I shall refer to the 
diatomic molecule as dihydrogen. We observe that even a simple molecule such as 
dihydrogen can exhibit a vast richness of properties. One such property is its nuclear 
spin isomerism [162] and forms a central topic of this chapter. 

The spin isomers are a result of the symmetrization postulate and the Pauli principle 
llb'3| 11571 H58| I164j . which impose certain symmetry requirements on a composite 
system comprising any number of identical particles: under the exchange of any two 
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particles, the total wavef unction must be either symmetric (even) or antisymmetric 
(odd). Furthermore, for particles, which are called fermions, the wavefunction must 
be odd under the exchange operation, also called the parity operation V p?T] . 

Suppose we have two fermions in distinct positions; I label their positions as A and 
B and assume their wavefunctions are x an d <f>- The overall wavefunction could then 
be written as if) = Xa4>b, and if we interchange their labels the resulting wavefunction 
becomes if)' = V{if>) = 4>aXb- Now, the Pauli principle simply requires the wavefunc- 
tion to be odd, i.e., if)' = —if). There is another interesting consequence of the odd 
symmetry of the fermionic wavefunction. If the fermions have the same position, say 
A, then 4>aXa = —Xa4>a implies a vanishing wavefunction, if)' = —if) = 0: two fermions 
with identical wavefunctions cannot occupy the same position, they are excluded from 
being in the same state — this latter fact being an expression of the Pauli Exclusion 
Principle. 

For dihydrogen, containing two protons, which are of course fermions, the overall 
wavefunction iptotai must be odd for the same reason. This wavefunction has contribu- 
tions from several degrees of freedom, all of which we assume act independently. These 
contributions come from the electronic states, translational, vibrational and rotational 
molecular motions and nuclear spin degrees of freedom, 



In keeping with the Pauli principle, I analyze the parities of the five independent con- 
tributions, and permit only those combinations for which the overall parity comes out 
to be odd. For the dihydrogen molecule the ground state electronic function if> e is even 
|164j . Exchanging the protons in the diatomic molecule has no effect on the transla- 
tional motion, which is simply a centre of mass motion, and so the translational part 
iptrans is also even. Likewise, for its vibrational motion, the molecule can be modelled as 
a pair of identical masses connected through a spring; interchanging the masses leaves 
the system unchanged, and so if) v yo is also even. The first three contributions in (|2.1j) 
being even, the parity of iptotai is therefore determined by the combined parity of the 
last two terms, if) ro t and if) ns : even rotational terms must combine with odd nuclear spin 
terms and odd rotational terms must combine with even spin terms. The reasoning 
is simple: odd functions multiplied by even functions result in an odd product. This 
explains the exact correlation between the molecular rotation and nuclear spin degrees 
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of freedom. The natural question to ask next is how the individual wavefunctions ijj rot 
and ip ns behave under the exchange operation. 

The rotational wavefunction of the molecule is a representation |3T] of the rotational 
state vector in the position space 

<MM)= (9,<f>\J,mj), (2.2) 

where 9 and <fi are angles in space and J is the angular momentum quantum number 
of the molecule, mj being its projection along the quantization axis. The function 



?/Voi(#, 0) transforms under rotations like the spherical harmonic function Y r 



J 

mj 



and the parity operator rotates this function according to 

V(Aot(9,<f>)) oc Y^n-e^ + n) = (-1)^(0,0) cx (-l) J Aot(9, 0), (2.3) 

which clearly shows that the parity is controlled by the quantum number J through 
the factor (— 1) J : an even J (= 0,2,4,...) results in symmetric rotational functions 
and odd J (= 1, 3, 5, . . .) gives rise to antisymmetric functions. 

I now turn to the parities of the nuclear spin wavefunctions. The kets corresponding 
to i\) ns for two spins can assume four forms, 

|00) = IT,), (2.4) 
|11) = | Tlx), (2.5) 
1 (|01) + |10)) = |T ) and (2.6) 



-4(101) - |10» = |5 >. (2.7) 



y/2 
1 

71 

The first three of these states (j2.4H2.6j) are characterized by the spin quantum number 
5=1, and are referred to as the spin triplet states and are denoted in general as |T su b). 
The subscripts 1,-1 and correspond to the projection of the triplets along the 
spin quantization axis. The state \So), given in ([2.7)1 . is characterized by S = and 
is called the spin singlet. (Finding these nuclear spin vectors for a composite two- 
spin system is a straightforward application of the problem of addition of two angular 
momenta [HI])- The proton is a spin-half particle (a fermion) and so for dihydrogen 
Si = S*2 = 1/2, implying that the total angular momentum quantum number S can 
be 1 or 0. The projections for 5=1, are clearly ra\ = {1,-1,0} and for S = 
the only possible projection is m = {0}. Corresponding to the four distinct {S,m s } 
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Table 2.1: Correlation between molecular rotational and nuclear spin wavefunctions. 
Para and ortho spin manifolds of H 2 . 

pairs, four nuclear spin states are possible, which are exactly the states (|2.4j) - (|2.7|) . 
These kets also completely span the four dimensional Hilbert space of the two nuclear 
spins in dihydrogen and a basis comprising these basis elements is generally called the 
symmetrical basis. 

For the purpose of establishing the correlation with the molecular rotational states, 
we are interested in the symmetry properties of these nuclear spin states. The triplets 
are symmetrical with respect to proton interchange whereas the singlet is antisymmet- 
rical. For example, relabelling the spins in the singlet \S ) would result in 

\So) = ^(|01) - |10» ^ i(|10) - |01» = -\S ). (2.8) 

In short, the nuclear spin triplets |Ti), |T_i) and |T ) have even spin wavefunctions 
and combine with odd rotational states. The singlet \Sq) has an odd spin wavefunction 
and must combine with even rotational states. This correlation is summarized in Table 

o 

Now suppose we can sort the dihydrogen molecules according to their values of J, 
keeping all odd J molecules in one cylinder and all even J molecules in another. As 
molecular rotation and nuclear spin are intertwined, the "rotational sieve" also acts 
as a "spin sieve", separating the singlet and triplet nuclear spin states. As a result, 
the odd J cylinder will be an equal mixture of the three kinds of triplet molecules, 
called ortho and the even J cylinder will solely comprise singlet molecules — their spin 
states are of just one kind and the molecules are termed para. Therefore if we can 
somehow partition the even and odd J rotational manifolds, then in principle one of 
these compartments will contain a pure quantum spin state. This pure quantum state 
is the initial state of our quantum computer and will be discussed in subsequent parts 
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of this thesis. 



2.2 Para-hydrogen from the perspective of statis- 
tical mechanics 

How do we in practice, achieve the spin sorting described above? Remembering that 
zero is an even number, we can cool the molecules to their rotational ground state 
J = and expect all nuclear spins to be in the singlet state. This approach effectively 
freezes the rotation of the molecules, locking them into the rotational ground state, 
and as a result, traps the nuclear spins in the singlet state. It is like manipulating the 
spin manifold through a more congenial arbiter, the molecular rotational manifold — as 
we now explore. 

The rotational energy levels are populated according to the Boltzmann equilibrium 
distribution jSU I163| I164j , the fractional population n j in each state being given as, 

nj= ^(gjg s ex V (-Ej/kT)), (2.9) 

Z being the relevant partition function, Ej the energy of the state with degeneracy 
9j9Si k the Boltzmann constant and T the absolute temperature. In the rigid rotor 
model |164t [L65J of dihydrogen, the energy corresponding to the angular momentum 
Pj is, 

E.;= ^J= J ( J + 1 )^ ( 2 - 10 ) 
where I is the moment of inertia of dihydrogen. Defining the rotational temperature 

(El, 

Br = ^, (2.11) 

()2.9|) can also be written as, 

nj= ^{gjg s exp(-J(J + l)6 r /T)). (2.12) 

Among gases, dihydrogen has the largest rotational temperature of about 85 K (calcu- 
lated from data for rotational constants presented in fT62] ) . Furthermore, the energy 
difference between any two levels J and J + 1, 

(Ej+i - Ej)/k = 9 r (J + 1)( J + 2) - 9 r J(J + 1) (2.13) 
= 29 r (J + l) (2.14) 
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is proportional to the rotational temperature and J + 1. For example, in dihydrogen 
the energy gap between the lowest (J = 0) and first excited (J = 1) rotational states 
corresponds to a temperature of 2 8 r = 170 K. It is therefore possible to cool to low 
temperatures, ensuring that most molecules are in the J = state. All other diatomic 
molecules have smaller rotational constants and the <-> 1 gap is smaller, requiring 
substantially lower temperatures to preferentially populate the J = state. For ex- 
ample, in the homonuclear molecule comprising two deuterium nuclei, 2 B.2, the lowest 
rotational states are separated by 29 r = 86 K; furthermore, as the deuterium nucleus 
is a boson (nuclear spin 1 = 0), the lowest rotational (J = 0) state would be populated 
by ori/jo-deuterium molecules. 

In the absence of external electric and magnetic fields, the Jth rotational state is 
(2 J + l)-fold degenerate (gj = 2 J + 1), and the allowed projections are given as 

mj = -J, -(J - 1), . . . , J - 1, J, (J > 0). (2.15) 

We have also seen that for even J, only one spin state (the singlet) is allowed and for 
odd J, three spin states (triplets) are allowed, and so the spin degeneracy gs in (j2.9J) 
is one or three depending on the parity of J. These degeneracies allow us to write the 
partition function for the dihydrogen molecule, 

Z= (2J+1) exp (-J(J + l)9 r /kT)+ 3 (2J+l)exp(-J(J + l)9 r /kT). 

.7=0,2,4,... J=l,3,5,... 

(2.16) 

2.2.1 Para- ortho ratios 

With these expressions, we can determine the relative ratio of the para and ortho 
molecules 

Npara = Eleven + 1) eX P (- J( J + l)9 r /kT) 

Northo SJ2 J=odd (2J+l)exp(-J(J + l)9 r /kTy l ' ' 

The fraction is a function of temperature and can be calculated numerically for different 

values of T, however analytical expressions can also be found under certain assump- 
tions. For example, in the high temperature limit T ^> 9 r , the spacing between the 
levels is small compared to kT and the discrete levels can be treated as a continuum 
of levels, replacing the discrete sums in (|2.16p - (|2.17|) with continuous integrals. For 
example, consider the sum 

oo 

^2(2J+l)exp(-J(J + l)9 r /T). (2.18) 
j=o 
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Carrying out the substitutions J(J + 1) = x 2 , 2 J + 1 = 2 x dx and 9 r /T = a, the sum 
can be approximated as 

POO 

/ 2xexp(-ax 2 )dx (2.19) 

Jx=0 

which is in fact a standard integral |166j with the solution 

1 T 

~= T - (2-20) 

Therefore, the sum over only the even or only the odd J states (in the high temperature 
limit) is, 

ii=il. (2.21) 
2 a 2 6 r ( J 

These integrals allow us to calculate the para:ortho molecular ratio in the high tem- 
perature limit. From (|2.17|) and the values of the integrals, we get, 



T»e r T/(29 r ] 



1/3, (2.22) 



N ortho 3T/(26 r ) 

showing that at temperatures above 28 r , dihydrogen is essentially a ~ 1 : 3 mixture of 
para and ortho molecules. 

At lower temperatures, when T < 9 r , the sums in (|2.16j) - (j2.17|) must be calculated 
explicitly. The series rapidly converges as only the lowest J states (such as J = 0, 1, 2) 
are populated and it is both permissible and customary to take only the first few terms 



if 1 + 5 exp (-60 P /T) + 9 exp (-2O0 r /T) + 



(2 23) 

Northo 3 V3exp(-20 r /T) + 7exp(-120 r /T) + llexp(-3O0 r /T) + ...y v ' ; 

This expression can be evaluated numerically to estimate the percentage of para- 
hydrogen in a para- ortho mixture as a function of temperature. Results are presented 
in Table l2~2*l and Figure l2~Tl showing that at the temperature of liquid N 2 (77 K) we 
have roughly an equal mixture of para and orino-hydrogen and at temperatures at or 
below 20 K, we obtain essentially pure para-hydrogen. 



2.3 Preparation of para- hydrogen 

Table 12.21 immediately suggests a method of preparing pure para-hydrogen: cooling 
the dihydrogen molecules to around 20 K and waiting for sufficient time, expecting the 
para- ortho mixture to equilibrate into the pure para form. However, this picture is not 
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Table 2.2: Percentage of para-hydrogen as a function of temperature. 
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Figure 2.1: Percentage of para-hydrogen as a function of temperature T. The plot 
also indicates the para-hydrogen temperatures employed in our experiment [SH] and 
the previous work by Hiibler et al. |152j . These experiments are described in Chapter 
E3 The dashed line is a guide to the eye, representing a para-hydrogen fraction of 1/4, 
the limit at high temperature. 
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accurate, to say the least. Ortho to para conversion is normally forbidden by angular 
momentum selection rules |162[ 11631 1164j , and so we need a slightly more sophisticated 
approach. For ortho to para conversion, J must change by ±1 and at the same time 
the nuclear spins must reorient from the 5 = 1 triplet state into the S = singlet. 
This reorientation can take place for example, in the presence of dipolar couplings with 
neighbouring ortho molecules |167j . however the intrinsic conversion is extremely slow 
|167j and takes place over a matter of days. To expedite the equilibration between 
the ortho and para manifolds, we can introduce an exogenous source of a magnetic 
field gradient |168j , breaking the nuclear symmetry of the ortho molecules, and for this 
purpose a catalyst with paramagnetic properties is normally used. 

The ortho <-> para conversion takes place on or near the surface of the catalyst, 
rapidly achieving equilibration and the resulting para: ortho ratio reflects the temper- 
ature of the surface, in accordance with ()2.23|) . If the surface temperature is about 
20 K, the ortho to para conversion reaches almost 100% completion. Removing the 
catalyst and warming the para-hydrogen to room temperature does not significantly 
alter the purity of para-hydrogen as the para molecules do not convert back into the 
ortho form, by virtue of the same rules that had prohibited their generation in the first 
place. It appears as if the molecules "remember" the surface temperature, the ratio 
of the isomers continues to reflect the catalyst surface temperature 1 . Upon warming, 
the molecules redistribute themselves in the even J rotational manifold according to 
Boltzmann rules; however their nuclear spin state does not change from the antisym- 
metrical singlet. Theoretically, we can store the para-hydrogen for indefinite periods 
of time; practically, for hours and whatever little conversion takes place is either at the 
surface of the storage cylinder or due to traces of impurities including residual ortho 
molecules. 

Several experimental setups are being used for the generation of para-hydrogen 
in laboratories around the world (for an overview of the methods, see j!50j ). Our 
collaborators in York have devised a method that uses liquid He as the refrigerant, 
lr This suggests that the para- ortho ratio can act as a reliable temperature probe. For example, 
studying this ratio for the isotopomers of water 169 in different comets, nuclear spin temperatures of 
w 25 K have been suggested [1701 1171] . significantly cooler than the observed rotational temperatures. 
Due to the strongly forbidden ortho^para conversion, we can glimpse the inter-stellar environment 
billions of years ago — approaching the lifetimes of comets! 
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cooling dihydrogen to temperatures as low as 7 K, allowing production of pure para- 
hydrogen. The para-ortho equilibration takes place on an activated charcoal catalyst. 
The H2 outlet is equipped with valves that allow evacuation via a vacuum pump, 
and sample filling via a Young's NMR tube adapter. Furthermore, the apparatus 
is equipped with a thermocouple, which allows the temperature of the H 2 gas to be 
controlled (our studies use temperatures in the range 18 K to 20 K), as well as a 
hydrogen pressure regulator on the inlet line and a pressure gauge on the outlet. The 
apparatus is left on permanently, except when the liquid helium reservoir needs to be 
re-filled, ensuring that high-purity para-hydrogen is always available. 

The pure para-hydrogen results in extremely low spin temperatures £Q, effectively 
in the mK range, while physically the gas can still be at room temperature. The spin 
temperature is now decoupled from the lattice temperature and from the perspective 
of the spins, the lattice is at infinite temperature, an approximation I shall use in 
modelling decoherence in our spin system (in Appendix ??). In other words, the spin 
entropy is very close to zero while the thermodynamic entropy of the gas remains very 
high. In effect, this is like placing a tiny spin refrigerator inside a molecular oven, 
which is hotter by five orders of magnitude. 

2.4 Chemical addition reactions with para-hydrogen 

Despite the exceptional purity, the para-hydrogen singlet state is spectroscopically 
useless. The S — state is the solitary inhabitant of the even J rotational manifold 
and no other spin states are available to which nuclear spectroscopic transitions can 
occur |156j . So para-hydrogen is NMR silent. Furthermore, dihydrogen is a symmetric 
molecule and the two protons are magnetically equivalent [3]. There is no way to 
distinguish between the nuclei and so, if we wish to exploit the very high purity of our 
para-hydrogen product, we must somehow break the symmetry and transfer the singlet 
state to another more interesting molecule, which allows for some greater variety. It 
is therefore clear that some kind of chemical reaction is needed. In fact, chemical 
applications have steered research in para-hydrogen, which has now been used for a 
number of years, in probing the structure of transition metal complexes and the study 
of reaction pathways involving hydrogenation of unsaturated organic molecules. (For 
example, see |15Uj and the references therein.) 
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Figure 2.2: Photolysis reaction scheme for hydrogenation of Ru(C0)3(dppe) to produce 
Ru(H) 2 (C0)2(dppe). The Ph represents an aromatic phenyl group. The C and H 
atoms in the bidentate side chain are not shown. 



Figure 12.21 shows the reaction scheme for a photolysis reaction |172t I173j involving 
the addition of para-hydrogen to a precursor molecule, I, Ru(C0) 3 (dppe) where (dppe) 
indicates 1, 2-bis(diphenylphosphino) ethane |174j . The molecule is stable and does not 
normally add dihydrogen, but a UV photon knocks off one carbonyl group to generate 
an unstable intermediate II, which immediately reacts with dissolved para-hydrogen to 
give the desired product III, which I refer to as the dihydride. The intermediate, II, is 
a transitory species and is known to react with para-hydrogen almost instantaneously; 
the reaction, in fact, takes place on the sub-microsecond timescale |175j . The dihydride, 
III, contains two inequivalent hydrogen nuclei, H/ and H5, the former being trans to a 
carbonyl group and the other, trans to a diphenylphosphino group. The two hydrogen 
nuclei form an AX spin system in the presence of a strong magnetic field. 

However, only breaking the symmetry in the para-hydrogen molecule is not enough. 
The addition of para-hydrogen to the precursor I, must be coherent, preserving the 
singlet spin state, both in form and magnitude, quality and quantity. This means that 
the pure |5*o) state must remain unblemished during the chemical reaction, without 
dephasing into a mixture of states or decohering into the maximally mixed state, which 
is always accompanied by a loss of valuable polarization. The pair of hydrogen nuclei 
in the dihydride, in other words, must fully inherit the pure singlet state from the para- 
hydrogen molecule. This is possible only if the hydrogenation is pairwise, meaning that 
both protons bonding with a transition metal centre in a single molecule, are derived 
from the same para-hydrogen molecule. This concern for preserving the singlet is a 
major concern in our experiments involving para-hydrogen and will be discussed in 
detail in the following sections. 
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Figure 2.3: Compound IV is 1, 2-bis(diphenylphosphino) ethane and V is a "mixed" 
phosphino-arsino compound 1, 2-bis(diphenylphosphinoarsino) ethane 

In addition to the "dppe" compound I, we have also used an analogue compound 
IV (see Figure l2~3~|) . 1, 2-bis(diphenylarsino) ethane, replacing the P atoms with As. 
Collaborators in York are also investigating other variants of I and IV, such as V, in 
which only one P atom is replaced with an As, or deuterating the compounds in the 
phenyl groups or the CH side chain, either fully or partially. In the context of para- 
hydrogen enhanced quantum information processing, the search for new compounds 
with better decoherence properties and more amenable spectral features continues, 
and chemistry continues to guide our way. 

The reaction depicted in Figure 12.21 is photo-induced, which means that the re- 
action only takes place in the presence of a UV trigger. A very short (12 ns) flash 
of UV radiation (308 nm) from a pulsed excimer laser |172j is used to initiate the 
photochemical reaction. The reaction is triggered by the laser flash and is complete 
within one /zs |175j . This allows precise control over the reaction, which can now be 
initiated and terminated at precise times. The laser effectively behaves as a switch, 
preparing pure singlet states in an AX system, as and when they are needed. Instead 
of using a single burst of radiation from a pulsed laser, we can also apply radiation for 
an extended period of time, using a continuous wave (CW) laser. In fact, we have per- 
formed experiments using both the pulsed and CW lasers, the former, for purposes of 
quantum computation and the latter, for purposes of demonstration and investigating 
the time dynamics of the spin system, especially over extended periods of time. The 
CW experiments, for example, can be useful in studying the decoherence properties of 
the system. 
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2.5 Operator description of the nuclear spin states 
in the para- hydrogen experiment 



Investigating the para-hydrogen enhanced systems from a quantum mechanical per- 
spective [160, 1176] . requires a description of the nuclear spin states in the language of 
density operators. The allowed spin states for the dihydrogen molecule are given in 
(|2.4H2.7j) . Committing a slight abuse of notation, I drop the kets and use the same 
symbols to now represent the corresponding density matrices, 

71= iTiXTil, r_!= iT-iXTUl, T = \T )(T \ and S =\S )(S Q \. (2.24) 

Writing these density operators (I call them states 2 ) in the computational basis, we 
obtain, 



71 ={{1,0,0,0}}, T_x= {{0,0,0,1}} 
/n n n n\ / 



Tn 



and So 




110 
110 

\o o o oy 

These states can also be written in product operator notation as 





1-10 
0-110 





(2.25) 



T ±1 =^ + ^(±I Z ±S Z + 2I Z S Z 



Tq — — + —(2 I X S X + 2 IySy — 2 I Z S Z ) 
I4 1 

Sq = — h — (— 2 TcS^ — 2 TyS^ — 2 Jj-S^ 



4 



- (ZQ X - I Z S Z ) and 
+ (-ZQ X - I Z S Z ) = 



I4 

4 



(2.26) 
(2.27) 

IS. 

(2.28) 

Pure para-hydrogen is the nuclear singlet state and would therefore be represented 
by S given in l|23Bj) . 

I4 



Ppara — Sq — I" ( ZQ X I Z S Z 



(2.29) 



and pure oriao-hydrogen, being an equal mixture of the triplet states has an operator 
form, 



2 The exact meaning of the word 'state' should be apparent from context. 
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Similarly a mixed state with a singlet fraction F will have the form, 



p F =Fp para + (1-F)p, 



ortho 



= ^ + (^-^)(ZQ x + I z S z ), (2.31) 
which reproduces the states in (j2.29H2.3l1j) by setting F — 1 and F = respectively. 



2.6 Ex situ hydrogenation 

There are two kinds of experiments involving addition of para-hydrogen to a precursor, 
depending on when and where the reaction takes place [160. Ilfijj . The first method 
involves hydrogenating the precursor outside the magnet and slowly transferring it 
into the strong field of the spectrometer j!77j : slow enough to ensure that the trans- 
fer is quantum mechanically adiabatic [HI], but fast enough to be rapid in compari- 
son with the Ti recovery towards the Boltzmann populations. This experiment was 
called ALTADENA (adiabatic longitudinal transport after dissociation engenders net 
alignment) by its inventors j!77j . 

In ALTADENA, the dihydride is formed outside the magnet, spins I and S acquire 
almost indistinguishable frequencies Dj ~ uj$ => (uji — u) s ) <C 2irJ, and therefore 
constitute a strongly coupled A 2 spin system. Only the So state will be populated and 
the triply degenerate triplets {T ±1 ,T } will be separated from So by the J coupling 
|176j 1178] . (Appendix ??,) as shown in Figure E3ta) . The dihydride is then adiabati- 
cally transferred into the intense field with the result that the system becomes weakly 
coupled, (ui — us) 2nJ, forming an AX system. The eigenstates in the weakly 
coupled system are the product states and out of these four states, only |01) and 1 10) 
have a non-zero overlap with the singlet, 

|4=(01|01 - 10)| 2 = |— i=<10|01 - 10) | 2 = -, (2.32) 



V2 V2 



whereas, 



|-^(00|01 - 10)| 2 = |^(H|01 - 10)| 2 = 0. (2.33) 

Due to the adiabatic condition, the system ends up in one of these states and the 
details are described in Appendix ??. The experiment is schematically depicted in 
Figure I23fb). As only one state (assume \aj3) = 1 01) ) is populated, two transitions 
are allowed, one for each spin, as is shown in the Figure. What the diagram does not 
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faithfully depict is the fact that if dephasing and decoherence are ignored, then the 
population difference between the populated and unpopulated levels will be of the order 
of unity; the spin order and state purity is conserved and the signal strengths will be 
very high as compared to the thermal signal intensities, shown somewhat generously, 
in part (d) of the same Figure. 

The ALTADENA experiment is related in spirit, to the remote detection experi- 
ment |179j . in which the magnetization is encoded in a low-field environment and then 
the sample is physically transported to a high-field magnet, conditions optimized for 
detection. However, such techniques which involve shuttling the sample back and forth 
are instrumentally demanding and for purposes of quantum information processing, we 
can achieve the same spin purity using an experimentally simpler approach. 

2.7 In situ hydrogenation 

The ALTADENA approach relies on slowly varying the Hamiltonian so that the spin 
system always remains in an eigenstate of the new Hamiltonian. What if the Hamil- 
tonian changes suddenly, for example, by switching on the magnetic field B instanta- 
neously? This is the basis of a second kind of experiment, called PASADENA (para- 
hydrogen and synthesis allow dramatic enhancement of nuclear alignment), and was 
discussed in the original prediction |15(j| of signal enhancements that are achievable 
using para-hydrogen. The effect can be explained using the "sudden" approximation 
|HU] in quantum mechanics. Consider a state p(t) and suppose that at time t = 0, the 
Hamiltonian suddenly changes from Hi to H 2 . The approximation tells us that the 
state p{t = 0~) does not change at the very moment the perturbation takes place, i.e., 

p(t = 0+) = p(t = (T). (2.34) 
Furthermore, for t > + , the initial state p(t = + ) evolves under the new Hamiltonian 

pit > + ) = exp {-iH 2 t) p(t = + ) exp (+iH 2 t), for t > + . (2.35) 

In zero-field conditions, the state of the two spin system is clearly pi = p para (as in 
I2.29j) and whatever the Hamiltonian suddenly changes to, the new state of the system 
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Figure 2.4: PASADENA and ALTADENA experiments: (a) is the energy level dia- 
gram for the dihydride product outside the magnetic field. Only So is populated and 
the triplets are triply degenerate. Diagram (b) includes the ALTADENA population 
diagram and spectrum for small 9, ()2.43j) . The diagram assumes that u>s > wi, re- 
sulting in only the 1 01) state being populated. All spectra are drawn assuming J > 
and follow the sign conventions described in [41] . Part (c) shows the population dia- 
gram and spectrum for the PASADENA experiment with selective excitation, ([2.47)1 . 
Finally, (d) shows the population diagram and spectrum for the thermal state, ([2.39)1 . 
also serving as a reference for labelling the coherences. In all spectra, frequency is 
positively increasing from left to right and the vertical heights of the spectra are not 
to scale. 
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at time t = + remains unchanged, 

P (t = 0+) = p(t = 0-) = p pa ra = S . (2.36) 

Therefore p par a = So is indeed an accurate initial state and can be directly used for 
determining the dynamics of our two qubit computer after the Hamiltonian switches 
from one form to another. 

One means to instantaneously change the Hamiltonian is to carry out the exper- 
iment inside a strong field: the dihydride product is formed inside the magnet, in- 
stantaneously changing the spin system from A2 to AX and therefore changing the 
Hamiltonian. This is the scheme we employed in our experiments. An in situ hydro- 
genation, can be achieved, quite conveniently, with our photolytic reaction, already 
illustrated in Figure 12.21 in which a laser is used to trigger the formation of the dihy- 
dride, at times of our choosing. Furthermore, the reaction takes place sufficiently fast, 
justifying the validity of the sudden approximation, which dictates that the two spins 
before and immediately after the in situ hydrogenation are in the state Sq. This results 
in both the 1 01) and 1 10) levels being populated, as shown in the energy diagram in 
Figure EOT c). from which we deduce the possibility of four transitions, resulting in the 
characteristic PASADENA spectrum, a pair of anti-phase doublets. The exact form 
of the spectrum, however, depends on the details of the experiment and the excita- 
tion pulse, both of which I shall discuss shortly. It must also be remembered that 
population diagrams [4, such as the one shown in Figure l2~4T c). may sometimes lead 
us to the wrong impression: for example, it is not evidently clear whether the pop- 
ulation distribution sketched out in Figure I2.4f c) suggests an incoherent mixture of 
states, 1/2(|01)(01| + |10)(10|), or the coherent state 1/2(|01) - |10>)(<01| - (10|) = S . 
What we know from the PASADENA density matrix calculations, corroborated by the 
spectra, is that it is indeed the latter. 

2.8 PHIP signal forms and enhancements 

I have indicated that in ideal cases both ALTADENA and PASADENA preserve the 
purity of the quantum state. The perfect spin order of the pair of X H nuclei in the 
dihydride product, results in very strong signals, when compared with the signals de- 
rived from the thermal state. These exceptional polarizations, whether obtainable from 
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ALTADENA or PASADENA are often captured in the term para-hydrogen induced po- 
larization or simply PHIP 160 . In fact, signal enhancements in NMR experiments 
involving hydrogen were observed even before the theoretical prediction of PHIP was 
made, and they were wrongly attributed [180J to CIDNP effects. Puzzlingly, the en- 
hancements increased when hydrogen was stored in liquid N 2 temperatures for longer 
periods before the reaction, and only later, researchers realized the origin of these in- 
creased intensities and that storing the dihydrogen in ultra-cold environments simply 
increased the mole fraction of para-hydrogen. (For a historical review of the PHIP 
effect, see |16Uj .) 

The form and intensity of the PHIP signal depends upon the details of the ex- 
periment and the flip angle of the detection pulse. Consider a state p subject to a 
detection pulse the state transforms to p' = P^pPj and consider the signal vector, 
(|1.23j) . from p', which I denote Sg(p, P4). This vector indicates the relative directions 
(emission or absorption) and relative strengths of spectral lines; but additional infor- 
mation, including reference spectra, are needed to associate individual vector terms 
with individual spin resonances. 

Thermal signal 

For a two spin homonuclear system, the density matrix p t h in the high temperature 
limit kT >> hu can be approximated 3 as 

Ptn= ^- + ®(I Z + S Z ). (2.37) 

At 295 K and for a proton frequency of 400 MHz, the Boltzmann factor, B, is 6.48 x 10~ 5 . 
A 90 y hard pulse converts p t h to the state 

p' th = h + ^ {Ix + Sx ), (2.38) 

corresponding to the signal, 

Sg(p th ,90„) = |{l, 1,1,1}. (2.39) 

3 An equally valid approximation is I4/4 — B/A(I Z + S z ) in accordance with (|1.49|) : but the sign of 
B is immaterial for our purposes. 
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ALTADENA enhancement 

The ideal ALTADENA experiment results in a state (e.g.), 

Palt = |01)(01|= h + l -(I z -S z -2I z S z ), (2.40) 

which can be excited by a 9 y detection pulse to yield the observable terms, 

^(sm9) (I x - S x ) - l - (sin^ cos 9) (2I Z S X + 2I X S Z ). (2.41) 

If 9 = 7r/2, the vector becomes, 

Sg(p a , t ,90 v )= ^{1, 1,-1,-1}, (2.42) 

showing that the spectrum comprises two in-phase doublets, one being in positive and 
the other in negative absorption. If the excitation angle 9 is small, sin# rs 9 and 
cos 9 « 1 and the signal vector can be approximated as, 

Sg(p ait , small 0) « {0, 0/2, -0/2, 0}, (2.43) 

showing that only the inner two transitions will be excited, in accordance with the 
ALTADENA spectrum predicted in Figure l2~4Tb) . Comparing (J2HSH) and (J232J), we 
can compute the relative enhancement of the signal intensities, 

1/4 2 

*7«w = -^g = -g, (2-44) 

which is ~ 30864 for a proton frequency of 400 MHz at 295 K, showing that the 
ALTADENA signals will be about four orders of magnitudes stronger than the ther- 
mal signal. This remarkable signal-to-noise enhancement motivated research in the 
detection of minor concentrations of dihydride isomers, sometimes appearing as inter- 
mediates, in catalytic hydrogenation reactions. The enhancement is itself a proof of 
the coherent addition of the pair of 1 H nuclei to the precursor: the hydrogens add 
pairwise, both nuclei coming from the same dihydrogen molecule. 

PASADENA enhancement 

The state after the PASADENA style para-hydrogenation is given by the singlet state 
(J2. 2812. 2912. 36)) and to compute the enhancements, we must first consider a suitable 
detection pulse. A simple hard pulse will not work here: the singlet is an isotropic 
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state, with a deviation component l/2(— 2I X S X — 2I y S y — 2I Z S Z ) = — I • S, which 
is symmetric in the I and S spin operators. Hard pulses, also have symmetric spin 
Hamiltonians of the form I a + S a and therefore commute with the isotropic state, 

[I ■ S, 6(I a + S a )] = 0, where a = x,y, z. (2.45) 

So, the singlet remains invariant and therefore undetectable under hard pulses, sug- 
gesting the need to apply an asymmetric RF Hamiltonian. Soft pulses such 6iI y + 02S y 
with 9\ 7^ 02, asymmetrically excite the spin system and could therefore be used for 
detection. In our experiments, we use a soft 90 I y pulse, selectively exciting only the I 
spin. We can deduce the observable terms using the product operator notation, 

—I ■ S = — (— 2 I X S X — 2 I y S y — 2 I Z S Z ) 
¥ 2 (2 IzS x 2 I y S y 2 I X S Z ) 
^ l -(2I z S x -2I x S z ), (2.46) 

suggesting that the spectrum comprises two opposite anti-phase doublets. This is also 
captured in the corresponding signal vector, 

Sg(p para ,90/ y )= 1,1,-1}, (2.47) 

which suggests seeing a spectrum of the form {— , +, +, — } or {+, — , — , +} depending 
on the identities of I and S spin resonances, where "+" and "— " refer to positive and 
negative absorption lines respectively. The enhancement for the PASADENA 90/^ 
experiment is identical to the ALTADENA 90 I y detection, i.e., 

1/4 2 , 

*- = Ws = F (2 ' 48) 

Similarly, if we have a mixed state of the form (|2.31|) . the signal vector and enhance- 
ment will be rescaled by a factor of (4F — l)/3, indicating that the enhancements are 
maximum for pure para-hydrogen. 

In short, if we trigger the reaction inside the magnetic field using a short laser pulse 
and detect the singlet immediately after the product formation using a selective 90° 
pulse, we should observe two anti-phase doublets. Subsequently, allowing the spins to 
relax to the thermal state and detecting with a hard 90° excitation pulse, we shall see 
two in-phase doublets, and the intensities of the thermal signal would be smaller by a 
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Figure 2.5: Different modifications of the PASADENA experiment: (a) instantaneous 
PASADENA involves detection right after a short laser flash from a pulsed laser, (b) 
the delayed version of PASADENA when the singlet is allowed to evolve for a duration 
r before being detected, (c) incoherent PASADENA, representing hydrogenation for 
an extended period of time r h , which can be achieved using our CW laser, (d) isotropic 
PASADENA when extended hydrogenation takes place in the presence of an isotropic 
mixing sequence, shown here by an outlined rectangle. In each of these diagrams, 
the solid black rectangle represents a detection pulse of some kind, the grey rectangle 
represents laser irradiation, ACQ represents acquisition, and periods marked by e are 
infinitesimal delays, shown only for illustrative clarity. 
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factor of 2/B = 2hu/kT. This factor is ~ 30864 for a proton frequency of 400 MHz at 
295 K. 



As a part of a research collaboration, we have investigated four kinds of experiments 
exploiting the PHIP effect. 

1. In the instantaneous PASADENA experiment, we detect the singlet immediately 
after it is formed. The laser shot and the detection pulse are controlled through 
a common spectrometer trigger, synchronizing the two events. Instantaneous 
PASADENA is schematically sketched in Figure l2~5T a). 

2. We also performed a delayed PASADENA: detecting the state at a time r > 
after the dihydride product is formed. This is depicted in Figure l2~5T rA 

3. Another modification is incoherent PASADENA, which involves hydrogenation 
over an extended period of time r^, shown in Figure l2~BT c). 

4. The extended hydrogenation can also take place in the presence of an isotropic 
mixing sequence (to be described in Section l2.12|) . I call this experiment isotropic 
PASADENA, and is shown in Figure l2~5Td). 

To understand these different versions of the PASADENA, we must first look at 
the evolution of the singlet in the presence of strong fields. 

2.9 Singlet dynamics in strong field 

The dihydride AX spin system will evolve in the presence of the AX Hamiltonian, 



where f2/ and are rotating frame angular frequencies of the two spins. Defining 
average and differential frequencies, 



H = Sl T I z + Vl s S z + irJ 2 I Z S. 



(2.49) 



A 



£li + £ls 
2 

2 



and 



(2.50) 



(2.51) 
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State p [p, I Z S Z ] [p, I z - S z ] [p, I z + S z ] 



ZQ X ^0 

DQ y ^0 

I Z S Z 

Table 2.3: Important commutation relations of different states p with operator terms. 

the Hamiltonian can be rewritten as, 

H = (S + A)I Z + (E - A)S 2 + nJ2I z S z 

H c Hnc 
■A. 



= £(/ z + ^)+7rJ2/^ + A(/,-5 z ) (2.52) 
= F c + F nc , (2.53) 

where, in the last step, I have decomposed the Hamiltonian into two parts, H c and 
H nc , the first commuting and the second, non-commuting with the isotropic singlet, 

[# c ,I-S]=0 and (2.54) 
[F„ C) I-S]/0. (2.55) 

These Hamiltonians also commute between themselves, 

[H c ,H nc ]=0, (2.56) 

allowing us to make the simplification, 

exp (-iHr) = exp {-i{H c + H nc )r) = exp (-iH nc r) exp (-iH c r). (2.57) 

Important commutation relations between several operator terms are summarized in 
Table ESI 

Considering only the deviation terms (—ZQ X — I Z S Z ), the singlet So evolves for time 
r according to the prescription, 

-ZQ X - I Z S Z exp (-iH nc r) exp (-iH c r)(-ZQ x - I Z S Z ) exp (iH c r) exp (iH nc r) 
= exp (-iH nc r)(-ZQ x - I Z S Z ) exp (iH nc r) 

(as H c commutes with the singlet, from (|2.54j) ) 
= exp (-zA(J 2 - S z )r)(-ZQ X - I Z S Z ) exp (zA(/ 2 - S z )r). (2.58) 
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The evolution can be simplified even further, as the —I Z S Z term commutes with A(I Z — 

S z ), 

[A(I Z - S z ) , -I Z S Z ] = 0, (2.59) 

and so the problem is reduced to considering the evolution of the ZQ X term under the 
Hamiltonian term A(I Z — S z ). The discussion leads us to two important observations. 
First, the singlet does not evolve under the weak coupling 2 I Z S Z Hamiltonian and 
second, the evolution of the singlet is completely determined by the evolution of the 
zero quantum term ZQ X = I X S X + I y S y , and that too, under only the term A(I Z — S z ), 
while the —I Z S Z term remains invariant, 

-ZQ X - I Z S Z {exp {-iAr{I z - S Z ))(-ZQ X ) exp (iAr(I z - S z ))} - I Z S Z . (2.60) 

I now consider the evolution of only ZQ X under only A(I Z — S z ). Suppose we place 
the transmitter (Tx) frequency in the exact centre of the / and S spin resonances; under 
these conditions, Q$ — an d the angular frequency A becomes (Qj — (— Qj))/2 = 
Qi, which for simplicity, I denote as Q. With this definition, the spin resonances will be 
2 Q or 2tt8 radians s _1 apart, where I have defined S as the difference between the spin 
frequencies in Hz. (Frequency separation is 2Q = 2tt5 rad s _1 .) Now predicting the 
evolution of —ZQ X under the Hamiltonian A(I Z — S z ) is straightforward if we consider 
the commutation relation, 

[ZQ X , ZQ y ] = r Iz ~ S *\ (2.61) 



and, defining 4 , 



enables us to write 



(L - S z ) 

ZQ Z = 1 \ (2.62) 



[ZQ X , ZQ y ] = iZQ z , (2.63) 



which is in complete one-to-one correspondence with [I x , I y ] = il z . I rewrite the Hamil- 
tonian A(I Z — S z ) as 2 A(I Z — S z )/2 = 2AZQ Z , and predict the dynamics of ZQ X , 

ZQ X ZQ X cos (2 At) + ZQ y sin (2 At). (2.64) 

The evolution is illustrated in Figure 12.61 and shows that the ZQ X term evolves in 

a subspace spanned solely by zero quantum vectors, {ZQ X , ZQ y , ZQ Z }, continuously 
4 This definition is found in |181j but it is not widely used in standard NMR literature. 
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(a) (b) 

Figure 2.6: Correspondence between the dynamics of (a) single spin operators and (b) 
two spin operators in the zero quantum subspace. The correspondence is I z «-> ZQ Z , 
I x <-» ZQ X , I y <-> ZQ y and u <-» 2A. 

traversing the "transverse" zero quantum plane with a precession frequency 2 A = 2tt5, 
and sinusoidally interconverting between ZQ X and ZQy. 

2.10 Delayed PASADENA 

I have already derived the signal obtainable from instantaneous PASADENA, (I2.47J1 . 
and I now predict the outcomes from the other experiments. In delayed PASADENA, 
the pure singlet evolves for a time r under the Hamiltonian (|2.49|) , in accordance with 

So ~ -ZQ X - I Z S Z ^ - ZQ X cos (2 Ar) - ZQ y sin (2 At) - I Z S Z (2.65) 

= — ZQ X cos (2 ti5t) — ZQ y sin (2 n5r) — I Z S Z ~ p(t), (2.66) 

provided we place the transmitter in the exact centre of the resonances. (As mentioned 
in Section fl.1.21 "~" denotes the deviation part of the matrix.) The second of these 
equations, ()2.66|) . re-parameterizes the frequency in terms of the absolute frequency 
difference 5. The equation also shows that in time r = 1/(45), ZQ X fully converts 
into ZQy and if we wait four times longer, i.e., for r = 1/5, ZQ X would go full circle 
and end up where it started from. This points towards the possibility of stroboscopic 
observation: if decoherence is neglected, Sq would appear static at times which are 

multiples of the precessional time period, k/5, for k — 0, 1, 2, Furthermore, ZQ y = 

I y S x — I x S y is antisymmetric in the I and S spin operators, and so we can detect the 
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state with a hard pulse. In Section EOl I shall describe how we can use properly timed 
delays, fully converting the symmetric ZQ X into antisymmetric ZQ y and hard pulses, 
to effectively achieve a selective 90 I y pulse. 

The state p(r) (as in (|2.66jl ). can be detected using either a selective 90 I y pulse or 
a hard 90° pulse. The selective pulse, results in the observable terms, 

p(r) 90Iv,obs - > I z S x cos (2 ndr) — I z S y sin (2 ty5t) — I X S ZJ (2.67) 

which compares with ()2.46|) by setting (27t<5t) = 0. The resultant signal will take up 
the form, 

Sg(p(r),90/ tf ) = ^{-I,l,-exp(2 7r5r),exp(2 7r5r)}, (2.68) 

indicating that the I spin remains in the anti-phase absorption mode, while the S lines 
contain dispersive components, the phases varying sinusoidally with r, although they 
always remain 180° apart within themselves. A hard 90. detection pulse applied to 
p(r) results in the sate, 

p(r) (-I Z S X + I X S Z ) sin (2 ttSt), (2.69) 

and the spectrum, 

Sg(p(r), 90.) = Sin(2 4 ^ T) {l, -1, -1,1}. (2.70) 

2.11 Incoherent PASADENA 

In incoherent PASADENA, the CW laser irradiates the precursor for an extended du- 
ration Th- singlets are formed at different times and consequently, evolve by different 
amounts. If the duration ^> 1/5, the average state will be pi nc = —I Z S Z as the 
terms in the ZQ plane will have dephased completely, as shown in Figure 12.71 Inco- 
herent PASADENA is, in fact, a time-distributed ensemble of coherent processes, the 
averaging in time acts as decoherence in disguise 182 . 

The state pi nc is incoherent, yet the magnetic ordering is very high. Traditional 
PHIP experiments usually generate this state, after photochemical or thermal hydro- 
genation is allowed to proceed for a protracted length of time. A selective 90 I y pulse 
converts it into the state, 

-I Z S Z ^ -I X S Z , (2.71) 
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"transverse" ZQ plane 
(a) 



singlets formed at different 
times incur different phases 

(b) 



(c) 



Figure 2.7: The incoherent PASADENA experiment. As time t increases from (a) 
through (c), new singlets are continuously being created and previously formed singlets 
begin to spread out in the ZQ plane. At the end of the irradiation period, t = r^, the 
magnetization vectors are evenly distributed in the ZQ plane. 



resulting in an anti-phase multiplet on the I spin while S remains completely "silent" . 



Sg(p mc ,90/ y )= ^{-1, 1,0,0}, (2.72) 

and therefore could be used to identify the I or S identities of the multiplets. A hard 
90 y pulse does not render the state observable, —I Z S Z I X S X ; however, a 9 y pulse 
with 9 ^ 7r/2 converts pi nc into the state, 

Pine ~ —I Z S Z — (I z cos 9 + I x sin 0) (S z cos 9 + S x sin 6) 

= — (I z S g cos 2 9 + I Z S X sin 9 cos 9 + I X S Z sin 9 cos 9 + I X S X sin 2 9) 
^lsm(29)(I z S x + I x S z ), (2.73) 

with the resulting signal vector, 

Sg(p mc , B y ) = l - sin (20){1, -1, 1, -1}. (2.74) 

The spectrum is of the form { — ,+,—,+} or {+,—,+,—}, comprising two anti-phase 
multiplets, and the intensities are maximized for the angle 9 = tt/A. The 45 y hard 
pulse yields the spectrum, 

Sg(p mc ,45„)= -1,1,-1}. (2.75) 

o 

and the enhancement over the thermal signal ()2.39j) is found out to be, 

Vine = ^ = (2.76) 
70 



which is ~ 15432 for a proton frequency of 400 MHz at 295 K, and is half the enhance- 
ments in the ALTADENA flUBD and instantaneous PASADENA (l2~4 r Hl) setups. 



2.12 Isotropic PASADENA 

Isotropic PASADENA is a variant of incoherent PASADENA: hydrogenation is allowed 
to proceed for a duration r^, but in the presence of an isotropic mixing jHE] sequence. 
The experiment uses the CW laser and is shown in Figure l2*3T d). The singlet will evolve 
under the Hamiltonian ()2.49|) . However, consider replacing H with a new Hamiltonian, 

#j,weak = 7tJ2I z S z , or (2.77) 

^strong = 7TJ2I-S. (2.78) 

The state Sq commutes with these coupling Hamiltonians, and is therefore, a stationary 
state, ideally remaining unchanged, even for extended lengths of time. Isotropic mixing 
sequences, built around inversion (180^) pulses [3*5] . in general, achieve the strong 
coupling Hamiltonian, ()2.78|) . and therefore, hydrogenation in the presence of an ideal 
isotropic mixing preserves the singlet |152j . 

The Hamiltonian, (J2.78|) . comprises only a J coupling and no Zeeman (or chemical 
shift) terms. In this case, Qj = Q,s = 0, satisfying the so-called Hartmann-Hahn match- 
ing EJ I181j condition, under which the polarizations are mixed (shared) between the 
spins, much like coupled pendulums continuously exchanging their kinetic energies, 
while keeping the total energy constant |183j . This energy- matched polarization trans- 
fer is the basis of the total correlation spectroscopy, TOCSY |184j experiment, which 
allows polarization to periodically "hop" among spins, even if they are not directly 
coupled. 

The underlying concept in modifying a Hamiltonian is very simple and can be 
understood using a familiar example: the homonuclear spin-echo sequence jSlIZI, shown 
in Figure 12.81 The sequence, shown in part (a) of the figure, sandwiches a 180^ pulse 
between two periods of free precession under the weak coupling Hamiltonian, ()2.49|) . 
If the two spin system is observed after the second r/2 period, then the evolution can 
always be described by an effective Hamiltonian ifj jWCa k = 27rJ I z S z , which involves 
only the scalar coupling term. Successive inversion pulses can be applied at intervals 
of r, which should be less than |1/J|, and if the observation is synchronized with the 
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(a) 



(b) 



Figure 2.8: Modifying the AX Hamiltonian ()2.49j) . to make it "look like" H JjWeak = 
2nJ I Z S Z . The pulse sequence (a) is the spin echo and (b) is its equivalent description, 
provided observation is stroboscopic. The black rectangle in each case, denotes a 180-r 
hard pulse. The upper line is the spin / and the lower is S. 

repetition rate, the effective Hamiltonian will be H Jweak . 

It is possible to "sculpt" more complicated Hamiltonians, using recurrent combi- 
nations of pulses and interspersed delays of free precession, or back-to-back pulses, 
and the mathematical theory of coherent averaging (also called average Hamiltonian 
theory) [2*1 1185j . remains a powerful tool for analyzing such sequences. For example, 
different kinds of pulse sequences can be designed to attain an effective Hamiltonian 
of the form, ()2.78j) . Most of these isotropic mixing sequences use composite pulses 

"33 US] ) instead of the simple inversion pulses, as depicted in Figure l2~Bl The com- 
posite pulses can compensate systematic errors and act over larger frequency ranges 
without overly demanding power requirements. Different possibilities of 'composite 
pulse mixing' 5 sequences exist, and they are thoroughly investigated in |181j . 

Composite pulse mixing sequences are generally "built up" from their constituent 
pulses by a more or less algorithmic procedure, which is outlined in Appendix ?? for the 
MLEV-16 sequence [186J. I numerically investigated the performance of different naive 
and sophisticated isotropic mixing sequences, in terms of their ability to preserve an 
ideal singlet state. These numerical results are also presented in Appendix ??. The next 
chapter describes our experimental results from isotropic PASADENA demonstrating 
the approximate preservation of the singlet. 

5 I use this non-standard phrase in comparison with the established phrase of 'composite pulse 
decoupling' 1551 ITST] . 
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Experiment 



excitation pulse spectral form predicted enhancement 



ALTADENA 


90, 


{+, 


+,-,-} 


2/B 


Instantaneous PASADENA 


90/, 


{- 


+,+,-} 


2/B 


Delayed PASADENA 


90, 


{+, 


- - +} 


2/i3sin (2n5r) 


Incoherent PASADENA 


45, 


{+, 


-,+,-} 


l/B 


Incoherent PASADENA 


90/, 


{- 


,+,0,0} 


2/B 



Table 2.4: Summary of the main results for ALTADENA and the various kinds of 
PASADENA style hydrogenations. 



2.13 Summary of PHIP signals 

The main results including suitable detection strategies, spectral forms and signal en- 
hancements for different PHIP based experiments are summarized in Table 12.41 Pre- 
dicted spectra for the different experiments are presented in Figure 12.91 
The next two chapters discuss our PHIP experiments. 

1. Our primary motivation in harnessing the PHIP effect was to prepare pure initial 
states that could be used for quantum information processing. We experimentally 
demonstrate the preparation of such states using a single flash of a pulsed laser, 
utilizing instantaneous PASADENA detection. The resulting two qubit system 
is almost pure and also lies above the entanglement threshold. 

2. With the isotropic PASADENA, we can demonstrate interesting dynamical char- 
acteristics of the singlet state, including decoherence. 

3. Finally, we implement two quantum algorithms (Deutsch and Grover) using pure 
initial states, demonstrating for the first time the use of pure states in liquid state 
NMR for quantum computation; the experiments are taken up in Chapter 0] . 
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Figure 2.9: Predicted spectra from (a) an ALTADENA experiment with a 90 y detection 
pulse, (b) instantaneous PASADENA with 90 I y , (c) delayed PASADENA with 90 x , 
(d) incoherent PASADENA with a 45 y hard pulse and (e) incoherent PASADENA with 
a 90 I y pulse. The I spin is assumed to be towards the left of the spectrum. All of 
these spectra are drawn to the same vertical scale; with (a), (b) and (e) of similar size 
and twice as big as (d). In plotting (c), I assumed sin(27r5r) = 0.6, resulting in a 
signal smaller than (a) by a factor 1/0.6. 
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Chapter 3 

Preparing pure initial states using 
para-hydrogen 

Every computation must start off in a well defined initial state. As seen in Chapter 
one of the major objections against NMR as a scalable quantum computing technology, 
is the problem of initialization. This chapter deals with our initialization experiment, 
which is an important first and successful attempt at realizing an essentially pure 
quantum state, also capable of performing universal quantum dynamics. The state 
also lies above the entanglement threshold. 

In our initialization experiment, pure para-hydrogen is prepared. The magnetic 
order is then coherently deposited into a molecule through a precisely controlled, fast 
chemical reaction. The migration of spin order from para-hydrogen to the molecule of 
interest must take place coherently and without appreciable loss, and at the end we 
must determine the density matrix and thus its purity and entanglement of formation. 
The preparation must also take place in a single experiment, without the need for 
multiple runs (as in temporal averaging schemes for pseudopure state preparation |87j). 
Our experiment demonstrably achieves all these requirements. 

Section l3~Tl describes a method of filtering the state for quantum state tomography 
in a single experiment and Section 13.21 describes the actual pulse sequence used for 
detection. The experiment itself, with all accompanying technical details is surveyed 
in Section 13.31 The acquired spectra are then processed and interpreted; the data 
processing and the final results are presented in Section l3~H Finally, Section 13*31 shows 
some results from the isotropic PASADENA experiment, utilizing the CW laser. A 
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very concise description of this experiment can also be found in our published work, 

Unsurprisingly, the para-hydrogen route to achieving pure initial states has also 
been attempted previously |152j . although the purity reached in that case was e ~ 
0.1, well below the entanglement threshold. Hiibler and co-workers used the molecule 
Ir(H)Cl(H)(CO)(PPli3)2 as their quantum computer, which they obtained from the 
addition of a 1 : 1 para-ortho mixture to their precursor molecule IrCl(CO)(PPli3)2 
(Vaska's catalyst). Our results show an almost nine- fold improvement in the purity, 
and we attribute this increase to two factors. First, we use almost pure para-hydrogen, 
prepared at 20 K and using Equation (|2.31|) with singlet fractions F = 1 and F = 1/2, 
gives us a three-fold enhancement in the magnetic order straight away (see Figure 
12. lj) . Second, the original experiment permitted hydrogenation for an extended period 
of time, resulting in signal loss due to the rotating frame spin-lattice relaxation, Tip 
|187j (in the presence of the isotropic mixing sequence j!86j); and we believe relaxation 
diminished their purity by another factor of three. Our experiment overcomes this 
second problem by letting the reaction proceed to completion almost instantaneously. 
Pure para-hydrogen gas and faster reaction times, therefore, act together in bringing 
our purity levels very close to one. 

3.1 One-shot tomography 

Full state tomography requires multiple copies of p and so we have come up with a so- 
called partial "twirl" 1 operation, which prepares a state having only two independent 
parameters, and is therefore characterizable in a single experiment. Our procedure 
characterizes the state obtained after the assumedly perfect partial twirl, and not the 
initial state prepared by the laser flash. In this sense, the tomography is restricted to 
only a constrained state, in which several independent parameters have been zeroed 
out. (As I will discuss later, our experiments indicate that the constrained state is 
almost indistinguishable from the initial state.) 

The ideal initialization experiment, which is instantaneous PASADENA, produces 
the singlet state, S . If this were the case, no tomography is needed, as the output 

would be known before hand. However, more realistically, impurity terms mix in with 
1 The twirl is the subject of Chapter |SJ 
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S - This could be due to several possible reasons: the small but finite reaction times, 
ZQ X ZQ y mixing under the AX Hamiltonian, the appearance of Ti and T_i terms 
due to spin-lattice relaxation and reaction with residual ori/io-hydrogen. The partial 
twirl takes an input state, pi n , with all these error terms, and simplifies it into a known 
form, p out , 

partial twirl , , 

Pin > Pout- (3.1) 

A key requirement of the partial twirl is that it leaves the singlet unchanged, both in 
form and fraction. Our partial twirl satisfies this requirement and outputs an So/T /T m 
mixture (see Section ll.4.6|) . a state of the form, 

Pout = a S + bT + c (Ti + T_i). (3.2) 

The coefficients a, b and c are the singlet and triplet fractions, i.e., a = (So) = F, 
b = (T ) and c = (Ti) = (T_i). For a valid density matrix, a + b + 2 c = 1, acting as a 
constraint on the {a, b, c} trio; as a result, only two coefficients are linearly independent 
and one-shot tomography is possible. An alternative representation of p ou t is in the 
product operator-multiple quantum basis, 

Pout = ^ +pZQ x + qI z S z , (3.3) 

which parameterizes the state in terms of p and q, the two independent parameters, 
which are directly related to the singlet and T triplet fractions, 

p — — a + b and (3.4) 

g = l-2o-26, (3.5) 

and conversely, 

o = i(l-2p- 9 ) (3.6) 

b = - A {l + 2p-q). (3.7) 

For pure Sq, a = I, b = c = and p = q = —1. Tomography translates into determining 
the values of these coefficients. The initialization experiment therefore involves a single 
laser flash for making the dihydride, filtering it to bring it into the desired form, and 
finally measuring the state achieved; the overall scheme is shown in Figure 13.11 
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G z2 















l< 1/5 >! k 1/5 

Figure 3.1: Scheme of the initialization experiment. The dark grey rectangle represents 
a single flash from the pulsed laser, the black rectangle indicates a 90 I y detection pulse. 
The arcs show strong gradient fields, sandwiching a light grey rectangle, representing 
a 90 x pulse. The gradients and the intervening pulse comprise the partial twirl, also 
called the filtration sequence. 



3.1.1 Implementing the partial twirl 

The partial twirl comprises a pair of strong z gradients (crushes), G z \ and G z2 , of 
unequal strengths, surrounding a 90 x pulse. I assume all pulse elements to be ideal. 
We can understand the working of the [G z i 90 x G z2 ] sequence by following the fate 
of the various terms in the density matrix as the sequence progresses. Defining the 
intermediate states p\ and p 2 , the states transforms as, 

Gzl 90 x G z 2 /o \ 

Pin ► Pi > P2 > Pout- (3.8) 

I expand the deviation part of pi n in the product operator-multiple quantum coherence 
basis, 

Pin ~ ei(l X + Iy + S X + Sy + I X S Z + IyS Z + I Z S X + I Z Sy) 

(3.9) 

+ e 2 (I z + S z ) + e 3 (DQ x + DQ y ) + e 4 (ZQ y ) + {l(ZQ x ) + m(I z S z )}. 

This is not a completely arbitrary state, as several terms are tied up with the same 
coefficient; however, the state is sufficiently general for our purposes. The terms with 
the e coefficients are the error terms and must therefore be prevented from reaching 
the output stage, p out , while ZQ X and I Z S Z should survive the filtration. 

The first gradient crushes all the single and double quantum terms, shown with 
coefficients e\ and e 3 , resulting in, 

Pin ^ Pi ~ e 2 (I z + S z ) + e 4 (ZQ y ) + {l(ZQ x ) + m(I z S z )}. (3.10) 
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The Zeeman (I z + S z ) and anti-symmetric ZQ y terms are transformed into single 
quantum coherences under the 90 x pulse and are crushed by the second gradient, 

+ 5*2) e2(—I y — S y ) — ^» decohere, and (3.11) 

e 4 (ZQ y ) = e 4 (I y S x - I x S y ) ^ e 4 (I z S x - I X S Z ) decohere. (3.12) 

Finally, for the ZQ X and I Z S Z terms, the [90 x G z2 ] sequence mixes the I and m coeffi- 
cients, 

ZQ X = I X S X + lySy ^ I X S X + I Z S Z = \{ZQ X + DQ X ) + I Z S Z (3.13) 

^ l -ZQ x + I z S z and (3.14) 

hS z I y Sy = X -{ZQ X - DQ X ) (3.15) 

^ ~ZQ X , (3.16) 

=^ 1{ZQ X ) + m(I z S z ) 1(1 + m)(ZQ x ) + l(I z S z ) ~ p out . (3.17) 

If we compare the above output, p out , with ()3.3j) . we deduce that p = {I + m)/2 and 
q = I, immediately suggesting that if I — m — — 1, then p — q — — 1, showing that the 
singlet component in p in is completely conserved, as we desired. I now consider the 
effects of finite duration gradient pulses. 

3.1.2 Finite duration gradient pulses 

It is not possible to implement zero length gradients, so we can instead think about 
selecting the proper rather than just the shortest timings. For this purpose, we can use 
the results in Section 1231 ZQ X evolves in the transverse zero quantum plane according 
to (j2.64|l . and therefore, setting the duration of each gradient to t g = (1/5), we ensure 
that ZQ X is practically "frozen" in time. Furthermore, the I Z S Z term does not precess 
at all, and neither of these terms evolves under the scalar coupling. For a ball-park 
figure of 5 = 500 Hz, t g — 1/5 — 2 ms, and the complete filtration sequence lasts for 
about 2(1/5) = 4 ms, so no appreciable signal loss due to relaxation takes place in this 
small amount of time. (For example, considering a T2 of 500 ms, the signal decreases 
by only 0.8%.) 

In our initialization experiment, we use two schemes: one without (Figure l2~H a)) 
and the other with (Figure I3~T|) the partial twirl. State tomography shows very similar 
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results for the two experiments, indicating that the error terms (with coefficients e in 
(|3.9Jl ) are very small, and that the laser flash generates a state which is close to a 
Werner singlet. 

We can also adjust the length of the gradient pulses to determine the I and S spin 
identities of the spectral multiplets. Halving the duration of each gradient from (1/5) 
to (1/2(5) prepares the longitudinal state I Z S Z , 

ZQx IzS z * ZQ X I Z S Z = IxS x + lySy IzS z 

> IxS x + IzS z ~ lySy — DQ X + I Z S Z 

G Z 2,tg = l/28 



I z Sz, (3.18) 

which is converted to I X S Z by a 90 I y detection pulse, resulting in the spectrum, 

Sg(I z S z ,90I y ) = -1,0,0}. (3.19) 

The I Z S Z state formed above is the same state, (of course with a negative sign,) as pi nc , 
(|2.71jh and the partial twirl with the gradient durations reduced to half, is therefore, a 
way of making the incoherent state from the coherent state. A selective pulse renders 
observable only one spin, which is unambiguously assigned the label I and the signal 
()3.19|) together with ()2.47|) can then be used to determine the positive and negative 
amplitude directions in the spectral outputs. 



3.2 The detection pulse 

We have devised a detection pulse sequence, determining the parameters p and q, all 
that is needed to completely reconstruct the density matrix, p out . 

In fact, the detection pulse is simply the selective 90 I y pulse already described in 
context of singlet detection in the PASADENA experiment ( Section 12. 8|) . The singlet 
(with p — q — — 1) results in the spectrum ()2.47|) while the general state, ()3.3|) . with 
arbitrary p and q, results in, 

p(I x S x + IySy) + q(I z S z ) ° bs > p(-I z S x ) + q{I x S z ). (3.20) 

The coefficients do not mix and are seen on separate spins, the p coefficient prepended 
to the S observable and the q coefficient to the I observable. The resulting signal 
vector is, 

Sg(p out ,90I y ) = ~{q,-q,-p,p}. (3.21) 
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Excitation pulse P signal vector Sg(p out , P) 



90/, 



45. 



90. 



e 



y 



'y 



y 



y 



1/4 {q, -q,-p,p} 
l/8(p-q) sin(2 0){-l,l,-l 
{0,0,0,0} 
l/8(p -<?){-!, 1,-1,1} 



1} 



Table 3.1: Signal vectors from application of some hard and selective pulses on the 
output of the filtration sequence p out as given in ()3.3|) . 

The signals resulting from applying different kinds of detection pulses are summarized 
in Table O 

Two comments about the detection pulse are not out of place here. 

1. The zero order phasing of the NMR spectrum is arbitrary. The y phase of the 
detection pulse results in absorption Lorentzians but using for example, a 90 I x 
pulse instead would result in the dispersive spectrum, 



which could be easily brought into the absorptive form. 

2. Similarly we could have also applied a 90 S y instead of a 90 I y pulse without 
loss of generality. However, in our experiments and mathematical analysis, we 
consistently followed the 90 I y description. 

3.2.1 Implementing the selective detection pulse 

There is a repertoire of methods available in the NMR literature [U 1188^ I189j to ap- 
ply selective pulses to one spin in a multi-spin system, the commonest being shaped 
pulses. In all our quantum information processing experiments, we have achieved spin 
selectivity using sequences based on the Jump-and- Return method [2H EDI- These 
kind of sequences, originally developed for water suppression, use only hard pulses and 
interleaving delays. 



Sg(p out ,90I x ) = -{iq,-iq,-ip,ip} = iSg(p out , 90 I y ), 



(3.22) 
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Here I present an explicit construction for designing the Jump-and-Return imple- 
mentation of our 90 I y selective pulse. The I selective pulse can be written as, 



90 L 



Ahly 

45 5,, 



45 I y 
45 S_ 



(3.23) 



the top and bottom rows comprising pulses on the / and 5* spins respectively, and time 
ordering in the pulse sequences being from left to right. Now, we know that a 9 y pulse 
can be decomposed 2 as a composite rotation [U E2 H2] , 



allowing us to rewrite ()3.23|) as 



90/,, 



By = 


90, 9 Z 90-*, 


clS, 






= | 45/, 


90 I x 


45 I z 


~ 1 45^ 


90 S x 


45 S^ z 



(3.24) 



90 

90 S- x 



(3.25) 



which can be expressed in more compact notation as, 

45 y 90, 45±, 90_„ 



(3.26) 



where the sequence now comprises only hard pulses and "contra-axial" (±z) pulses. 
This sequence can in fact be simplified [TH| further. Rearranging the terms in (J3.24|) . 
we are able to deduce the relation, 



y 90, = 90, 6 Z 



(3.27) 



Furthermore, we can displace the z rotation to the left or to the right, while changing 
the phase of a neighbouring pulse. For example, using the identity, 



Pz = P. 



z tX(/>+P 



(3.28) 



the sequence (|3.27|) becomes, 



e y 90, = 90, e z = e z 90 0+e = e z m b 



(3.29) 



This allows us to re-write ()3.26|) in the form, 



45 2 90 45 45 ±z 90. 



(3.30) 



2 This identity can be easily proven by explicit calculation. 
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a sequence comprising a non-selective z rotation, a pulse, a contra-axial ±z rotation and 
a final pulse; both pulses being hard. This sequence is a state independent prescription 
for implementing a 90 I y selective pulse. For the output of the filtration sequence p ou t 
()3.3j) . comprising only ZQ X and I Z S Z terms, a further simplification is possible, as both 
terms in p out commute with non-selective z rotations. In fact, the very purpose of the 
transformation (|3.29|) was to bring this z rotation to the beginning of the sequence; it 
has no effect on the state p out and could therefore be simply dropped altogether. The 
pulse sequence for fully characterizing the filter output can therefore be written as, 

90 45 45 ±z 9CU, (3.31) 

two hard pulses sandwiching a properly timed delay, in the spirit of the Jump-and- 
Return style selective pulses. The hard pulses are straightforward to implement, while 
the 45± 2 rotation is achievable by observing its propagator form, 

exp(-i~(I z -S,)). (3.32) 

From this, we immediately recognize that a suitable Hamiltonian to implement this 
function would be the non-commuting part H nc = A(I Z — S Z ) of the background Hamil- 
tonian H as described in ()2.53|) . This Hamiltonian is achieved by necessarily placing 
the transmitter in the exact middle of the I and S spin resonances and allowing the 
two spin system to evolve under the background Hamiltonian H ()2.49|) for a time t^, 

71 7T TT 1 , 

At *=4 =- <*= *z= m= (3 - 33) 

For a representative value of 5 = 500 Hz, this delay is 500 ps. 

The complete pulse sequence for the initialization experiment can therefore be writ- 
ten as, 

flash G zl) [ tg ] 90 x G x2 ,[t g ] 90 45 [td] 90-x acquire, (3.34) 

schematically represented in Figure l3~2T a). If the laser flash produces a hydride without 
any error terms in p in , we could rather leave out the partial twirl altogether; the 
resulting scheme is shown in part (b) of the figure. We have performed both kinds of 
experiments and observe no appreciable difference between the results. Interestingly, 
if we knew that p out is indeed the singlet with p = q = — 1, we could also leave out the 
first hard pulse in the tomography sequence (j3.31|) as the singlet commutes with hard 
pulses. However, for mathematical completeness, we retain this initial pulse. A few 
additional comments about the detection pulse sequence are the following. 
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hv 90 x 90 45 90. x hv 90 45 90. x 




ACQ 



(a) (b) 

Figure 3.2: The initialization experiment pulse sequences with (a) and without (b) the 
partial twirl. All rectangles are appropriately labelled. The transmitter frequency is 
placed in the middle of the resonances. 

1. Analysis shows that the inter- pulse delay, = 1/(4 5) maximizes the ZQ X to 
ZQ y conversion, acting as a tt/2 rotation in the transverse zero quantum plane. 

2. The sequence is most robust against pulse width errors when the phases of the 

pulses are 135° apart. The role of errors in this sequence is considered in Appendix 
?? 



3.3 Experimental system and methods 

In this section, I describe the experimental details of our instantaneous PASADENA 
experiment 3 . 

3 A11 experiments discussed in this and the next chapter, were performed by Dr. S. B. Duckett 
and Dr. D. Blazina at the University of York. This included synthesis and characterization of all 
compounds, manufacture of the para-hydrogen plant, building and adjusting the optics for efficient 
transfer of light to the sample, and finally implementing the pulse sequences on the spectrometer. 
Myself and my supervisor, at Oxford, were involved in pulse sequence design and analysis, data 
processing, state tomography, decoherence modelling, simulations and the theoretical design of the 
experiment. However, I visited York a couple of times to assist in the experiments and also carried 
out local experiments on simulated states. 
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III 

Figure 3.3: The product of the photolysis reaction, Ru(H) 2 (C0)2(dppe). The two 
qubits in the quantum computer are the protons H/ and H5. 

3.3.1 Chemical system 

The molecules used in our experiment have already been identified in Section 12.41 and 
the photolysis reaction underpinning pure state preparation is illustrated in Figure l2~2l 
The stable precursor, I, Ru(C0)3(dppe) is photolyzed using a laser flash, converting 
it into the transitory species, II, Ru(C0)2(dppe), which immediately reacts with dis- 
solved para-hydrogen producing Ru(H) 2 (C0) 2 (dppe), III fFigure l3~3*j) . the molecule of 
interest. The reaction takes place inside a magnet and is an example of instantaneous 
PASADENA. 

The dihydride III is a well-characterized system [174j and no secondary isomers 
are known to be produced during or after the photolysis. The reverse reaction, III — > I 
is extremely slow at the operating temperature of 295 K. (The backward reaction is 
catalyzed by higher temperatures, above 340 K, and by the presence of C0.) 

The precursor I, was dissolved in deuterated benzene (CqDq) and placed in a 5 mm 
(diameter) NMR tube; from which dissolved gases were subsequently removed using 
freeze-pump-thaw cycles, corresponding to, respectively, freezing with liquid N 2 , de- 
gassing and finally melting the sample. The tube was covered in Al foil to protect 
against premature photolysis by ambient light. Para-hydrogen prepared at 20 K was 
introduced into the tube, at a pressure of about 3 bar. Warming and subsequent shak- 
ing ensured that the p-H 2 gas dissolved in the solvent. All NMR studies were carried 
out with samples of approximately 1 mM concentration on a Bruker DMX-400 spec- 
trometer with 1 H at 400.1 and 31 P at 161.9 MHz, respectively. NMR properties of the 
system and several other experimental parameters are summarized in Table 13.21 
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System NMR parameters 
Ru(H) 2 (C0) 2 (dppe) 



resonance frequencies —6.32 and —7.55 ppm 

5 492 Hz (on the 400 MHz spectrometer) 

: H Ti (single quantum) 1.7 s 

*H T 2 (single quantum) 0.58 s 

2 J HH 4.6 Hz 



Experimental parameters 



Proton frequency 400 MHz 

System temperature (T) 295 K 

t n / 2 8.25 /is 
Number of laser flashes (F) for control * 1000 

Inter-flash delay (tf) * 7 s 

Number of scans (S) * 3072 

Inter-scan delay (t s ) * 20 s 

Active volume fraction (V» * 12.5/34 = 0.368 

Spectral width ps 30 ppm 

Receiver gain 128 



Table 3.2: Important system and experimental parameters involved in the initialization 
experiment. Parameters marked with a "*" are described in the text. 
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(b) 

Figure 3.4: One dimensional image of the (a) NMR tube filled with water and (b) of 
the dihydride after the laser flash; showing that it is indeed formed within the coil 
region. The "hump" in the active region is mainly due to non-deuterated protons in 
the benzene solvent. 



3.3.2 Optical assembly 

The spectrometer was fitted with a special probe equipped for in situ photolysis |173j . 
also designed by colleagues at York. For the single flash experiment, an MSX-250 pulsed 
XeCl excimer laser radiating 308 nm UV light was employed. The flash lasted for 12 ns 
and the flash-to-flash variability was known to be less than 1 ns. The specification sheet 
of the laser quotes a beamwidth of 12 x 6 mm and a beam divergence of 2 x 2 mrad; and 
the power per pulse was measured to be 32 mJ. The single pulse of UV light was then 
directed into the probe of the spectrometer using a carefully aligned assembly of quartz 
prisms, an arrangement resembling a periscope. The UV light was directed to fall in the 
active region of the NMR tube, lying directly below the RF coils; this is made possible 
by precise engineering of the probe head. The laser is triggered by a control signal, 
initiated from within the pulse programme. There is also a small (~ 10 /xs) trigger-flash 
delay and the tomography pulse or filtration sequence can "wait" to take into account 
this small time interval. The dihydride is formed within the coil region because of 
the optical positioning; this is also confirmed by a one dimensional image [HH] of the 
dihydride taken after the laser flash, showing the relative positions of the hydride 
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and the active region and is shown in Figure 13.41 The experiment was performed 
in a darkened room, reducing the possibility of background photodissociation of the 
compound, ensuring that all dihydride formation is initiated by the UV laser. 

For extended durations of irradiation, a Kimmon IK3202R-D 325 nm He-Cd 27 mW 
CW UV laser was used; the optical assembly was left unchanged. To enable a pre- 
determined amount of laser light to reach the sample, the CW laser was employed in 
conjunction with a Vincent Associates Uniblitz VMM-D1 shutter, which was opened 
and closed from within the logic of the pulse programme. The response time of the 
shutter was 1.5 ms; and since our CW experiments involved minimum irradiation times 
of 60 ms, errors due to shutter timing were negligible. 

3.3.3 Dihydride spectra 

The NMR spectra for different pulse sequences, in conjunction with a single flash 
of the pulsed laser, are shown in Figure 13.51 The GARP broadband decoupling se- 
quence (190J was applied throughout signal acquisition to remove couplings to 31 P 
nuclei. For the initialization experiment, the NMR tube must not contain any resid- 
ual Ru(H) 2 (C0) 2 (dppe) before the flash is triggered. This is verified by applying a 
simple 90 y hard pulse prior to irradiation, no signal is seen, suggesting that no dihy- 
dride exists before the flash or its concentration is too minute to be detectable. The 
resulting spectrum, which is essentially noise, is shown in part (a) of the figure. Figure 
I3.5f b) is the result of performing instantaneous PASADENA experiment depicted in 
Figure l3~27 b). a selective Jump-and- Return sequence applied soon after the flash. The 
spectrum shows the expected pair of anti-phase doublets, in accordance with ()2.47|) . 
Part (c) of the figure, shows the spectrum resulting from the sequence (|3.34j) . but this 
time also includes the filtration sequence before the Jump-and-Return. The result is 
very similar, and in fact apparently indistinguishable from (b), suggesting that the 
para-hydrogen addition is clean, and that error terms mixing in with the singlet are 
very small. Modifying the filtration sequence by reducing the duration of each gradi- 
ent to 1/(2 5), results in only the / spin being observed. Part (d) of Figure l3~B1 shows 
exactly this, authenticating the right multiplet as belonging to I. With the identity of 
the spins now revealed, we can translate the spectral pattern in Figure 133T b). (read 
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Figure 3.5: NMR spectra from the single flash experiment. The spectrum in (a) is the 
output from a 90 y hard pulse before the laser flash; (b) is the spectrum from the pulse 
sequence (I3.31j) . i.e., without the filtration; (c) is acquired after the sequence (|3.34|) . 
inclusive of the filtration; while in (d), the length of each gradient is halved, exciting 
only the / spin; and finally, (e) is the spectrum obtained from the sequence (13. 31 J) , but 
without the initial 9O45 hard pulse. 

leftwards as down, up, up, down), into the signal vector { — , +, +, — } 4 , where the left- 
most entries in the vector correspond to spin / and the rightmost to S. For example, 
with this convention the spectrum (d) translates to {+, —,0,0}. In short, the pair of 
lines appearing to the right in our displayed spectra belong to spin / and correspond 
to the two leftmost positions in the signal vector. Finally, part (e) of the figure shows 
the spectrum obtained from the pulse sequence [tj, 90_ x ], dropping out the first hard 
pulse from ()3.31|) and the result is apparently indistinguishable from (b). Hard pulses 
are irrelevant for the singlet, indicating that the state prepared after the laser flash is 
indeed close to being a perfect Werner singlet. Our tomography results in Section EOl 
will in fact, state the exact component-wise breakup of the quantum state. 

FigureHnHshows spectra from the delayed PASADENA experiments (Figure lTBT b)). 
in which the selective detection pulse is postponed for a variable delay r after the laser 
flash. The results are in complete agreement with (|2.68|) . with the / multiplet remaining 
more or less unchanged (neglecting relaxation), and the phases of the component lines 
in the S multiplet varying sinusoidally, with the expected frequency of S. 
4 I assume '+' is upwards and '— ' is downwards. 
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Figure 3.6: Spectra from delayed PASADENA, in which the selective detection pulse 
sequence is initiated with a variable delay r after the laser trigger. The plotted spectra 
were obtained by systematically varying r under computer control: increasing in steps 
of nearly 1/(45) ~ 510 /xsec, from 3 fis in (a) to 2043 fis in (e). Multiplet I remains 
unchanged (ignoring relaxation), while the components of the S multiplet change phase 
in steps of 7r/2; going from (a) absorption {—,+,+,—}, to (b) dispersion {— , +, i, — i}, 
(c) negative absorption {— , +, — , +} to (d) negative dispersion {— , +, and back 

to (e) absorption {—,+,+,—} . This cyclic variation of the phase occurs at a frequency 
6. 
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3.3.4 Control experiment 

The experimental spectra (Figure E3J) are in almost perfect agreement with the theo- 
retical predictions. However, we are interested in initializing the two qubit system in a 
pure state, and to find the exact purity, it is essential to know not only the form, but 
also the intensity of the signal. As absolute intensities are hard to come by in NMR, 
we need a reference spectrum with a known intensity. The thermal signal, specified 
in (|2.39|) . is a readily available reference with an intensity deducible from first princi- 
ples. Our control experiment provides this intensity reference and thereby calibrates 
the enhanced hyper-polarized signal strengths. 

One straightforward means to implementing the control is allowing the highly po- 
larized dihydride to relax to the thermal equilibrium state (|2.37j) and applying a hard 
90^ pulse, followed by acquisition. Indeed this is a completely valid strategy, but the 
problem is that the amount of dihydride, Ru(H)2(C0)2(dppe), formed from a single 
laser flash is much too small and, after relaxation, the thermal signal is too weak to be 
detected. For example, assume we make a pure state in our initialization experiment, 
then if the thermal signal from a single laser flash is to be detected over and above 
the noise floor, the signal-to-noise ratio for the pure-state spectrum should be of the 
order of 30000: this is something clearly not possible. Our spectra in Figure 1331 show 
a signal-to-noise ratio of only about 40. 

One way around this limitation is to increase the concentration of the dihydride by 
applying more flashes. The initialization experiment with F flashes 5 is schematically 
shown in Figure 13*771 In practice, we applied 1000 laser flashes in total, thereby increas- 
ing the concentration of the dihydride a thousand-fold. Even then, the signal-to-noise 
of the thermal spectrum was not sufficient, and it was necessary to perform signal 
averaging We performed 3072 scans before a thermal spectrum with a reason- 

able signal-to-noise ratio was obtained. The control and enhanced spectra are shown 
side-by-side in Figure 13.81 The resulting thermal signal is a sum of S scans and for 
meaningful comparison with the para-hydrogen enhanced spectrum from a single flash, 
must be divided by S x F. We can then compare the para-hydrogen with the properly 

rescaled thermal spectrum and determine the enhancement and hence the purity. This 
5 The notation F should not be confused with the identical notation I have used for the singlet 
fraction. 
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Figure 3.7: Schematic outline of the calibration experiment. The first laser flash, 
shown by a grey rectangle, produces Ru(H) 2 (C0)2(dppe) which is detected and the 
signal acquired (A). Subsequently we apply F laser flashes producing F times as 
much dihydride. Adjacent flashes are separated in time by tf. The thermal signal was 
measured by signal averaging S times (A\ through A s ), and the averaged spectrum 
obtained after signal averaging, is in fact, the sum of the S transients. Two successive 
scans are separated by a delay of t s . 



naive approach, however, yields an enhancement of about 77000, well above the the- 
oretical maximum. The apparent discrepancy is addressed and resolved by an active 
volume argument and is the subject of the next subsection. 

In the preceding discussion, we have assumed that F flashes produce F times as 
much Ru(H) 2 (C0)2(dppe) than the amount produced by a single shot, or in other 
words, the precursor I is not depleted. Our collaborators, in fact, performed a 31 P 
quantitative experiment [35] to examine the extent of depletion, which turns out to 
be very small. The underlying concept of the quantitative experiment and results are 
outlined in Section 13.3.61 

An inter-flash delay tf of 7 s is sufficiently long to allow different precursor molecules, 
which are subject to convective and diffusive motions, to cross the path of the laser 
beam, and be irradiated by successive laser flashes. This ensures that there is, more 
or less, a linear relationship between the amount of Ru(H) 2 (C0) 2 (dppe) formed and 
the number of flashes |153j . As a matter of principle, the selected tf should be large as 
compared to the characteristic timescales for diffusion and convection. Furthermore, 
successive scans are separated by about 20 s, which is greater than 5Ti « 8.5 s and so 
saturation effects [33] can also be safely ignored. With these values of tf and t s , the 
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Figure 3.8: The para-hydrogen enhanced and calibration spectra of 
Ru(H) 2 (C0) 2 (dppe). The enhanced spectrum (a) is a single scan after a single 
laser flash, and is identical to Figure 133T b) . while the calibration spectrum (b) is the 
sum of 3072 scans after 1000 laser flashes. The calibration spectrum has been divided 
by 3072000 and then multiplied by 30864 (the theoretical maximum enhancement) 
before plotting, so that the two spectra should show the same intensity. In fact the 
para-hydrogen enhanced spectrum (a) is even more intense than naively predicted. 

total time required for the experiment is about 19 hours. 
3.3.5 Active volume correction 

Figure 13.51 shows the spectrum from a single laser flash alongside the calibration spec- 
trum, which has been multiplied by 30864/ (S x F) before plotting. For an almost pure 
state the two spectra should be roughly the same size, but the para-hydrogen enhanced 
spectrum is significantly bigger than the thermal spectrum, resulting in an apparent 
enhancement about 2.5 times larger than the maximum possible. This discrepancy is 
explained by considering the fact that the NMR probe is not sensitive to the entire 
sample, but only to an "active" region, lying inside the RF coil (refer to Figure l33j) . 
Only the sample that lies within the active volume can be detected, the rest remains 
undetected, and therefore lies in what I call the "passive" region. 

A single laser flash produces Ru(H)2(C0)2(dppe) lying wholly within the active 
region and therefore, all of it is detected and contributes to the enhanced signal. This 
is not true for the calibration signal. Applying the 1000 laser flashes takes about 2 
hours during which the dihydride distributes itself throughout the entire volume of the 
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sample by diffusion and convection. As a result, the calibration spectrum does not 
record all of the dihydride present in the NMR tube, but only registers the fraction 
lying inside the active region; molecules in the passive region yield no spectrum. As a 
result, the calibration signal is attenuated by a factor Vf, the active volume fraction, 
corresponding to the ratio of the active to the total volumes. The true enhancement 
is therefore smaller than it appears, the correction factor being Vf. As expected the 
directly measured enhancement (that is, before correction) shows a linear dependence 
on the total sample volume (data is shown in our submitted work |153j ). confirming 
our active volume hypothesis. 

The active volume fraction can be estimated using simple geometrical considera- 
tions. Considering a ratio of lengths as a good substitute for a ratio of volumes, Vf 
can be estimated from the length of the RF coil and the depth of the NMR sample. 
In the experiment described, these lengths were measured 6 to be 12.5 mm and 34 mm, 
resulting in an approximate fraction Vf of 12.5/34. This deduction of the Vf fraction, 
however, is prone to error, as it is based on an "all or none" assumption of the probe 
sensitivity, that is, all of the dihydride inside the active region and none in the passive 
region is detectable. But as Figure I3.4f a) shows, the sensitivity drops off smoothly 
near the edges of the RF coil, rather than in a sharp step-like manner. As a result, 
the sensitivity is not uniform even within the active region, and our estimated Vf re- 
mains, at best, only an estimate 7 and addressing this should be a priority in future 
experiments. 

3.3.6 Depletion of the precursor 

The control experiment assumes that all flashes behave identically: having identical 
photon fluxes, quantum yields and more importantly, that the precursor (I) concentra- 
tion is constant at all times during the experiment. This section considers the effects 
of depletion of the precursor molecules, ignoring any variability in the light pulses. 

Under this assumption each flash delivers a fixed number of photons to the sample. 
Suppose this number is R. Each of the R photons can potentially affect the conversion 
of a precursor molecule to the dihydride, however the number of molecules converted 
6 The sample is long because this improves the homogeneity of the magnetic field by moving the 

susceptibility boundary away from the detection region. 

7 My rough guesstimate for the error in the active volume fraction would be 10-30%. 
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depends on the concentration of the precursor, [C], existing just before the flash. A 
simple first order model for the conversion can be written as, 




where A is a rate constant independent of the concentration. What this constant 
signifies is that although the absolute number of photons converted will fall with each 
successive flash, the fraction of the molecules converted per flash remains constant. For 
example, the thousandth flash would convert a smaller number of precursor molecules 
to the dihydride than the first pulse, but in each case the fractional conversion would 
remain the same. Now suppose that each flash converts a constant 0.01% of the precur- 
sor. If depletion was not an issue, then 1000 flashes would convert 1000 x 0.01%= 10% 
of the precursor, but taking depletion into account, we predict an overall conversion 



which is 9.52%, smaller than the 10% as expected. Thus the no depletion assumption 
underestimates the thermal signal and overestimates the purity of the para-hydrogen 
state by a factor of 10/9.52 = 1.05, which is our depletion correction factor. In the 
hypothetical limit of an infinite precursor concentration, the rate of conversion d[C]/dt 
becomes more or less independent of [C] and the conversion acts as a pseudo zeroth 
order reaction in which the absolute number of molecules converted per flash appears 
to be constant. This is because the initial concentration of the precursor is so high 
and the conversion per flash so low, that [C] does not substantially decrease on the 
application of successive flashes. Mathematically, this implies that all three parameters 
in the R.H.S. of ()3.35j) act like constants, resulting in a linear decrease in the precur- 
sor concentration (or linear increase in the dihydride concentration). These ideas are 
illustrated in Figure 13.91 

In fact, we carried out a simple experimental calculation to determine the extent 
of precursor depletion. The two phosphorus nuclei acquire different chemical shifts in 
I and III, |174j . resulting in distinct 31 P signals for these compounds. This fact helps 
in performing a quantitative experiment 8 , determing the relative concentrations of the 
8 The experiment involves direct 31 P detection; uses 10000 scans (the large number of scans are 
needed because of the low concentrations of the precursor), an inter-scan delay of 20 s and proton- 
decoupling turned on only during acquisition, preventing build-up of artificial intensities due to the 
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(a) (b) 

Figure 3.9: Depletion of the precursor: in (a), the solid line shows the overall fraction 
of molecules converted, (1 — (1 function of F. The dotted line is a plot 

of Fx, the total fraction of molecules converted for a hypothetical pseudo-zero order 
conversion and coincides with the depletion curve for small values of Fx. Part (b) 
shows the depletion correction factor, Fx/(1 — (1 x) ) ? clS db function of F. The 
fraction converted per flash x is assumed to be 0.01%. 
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two species after F laser flashes, and therefore, the extent of conversion of I to III. One 
such experiment (see footnote) showed that after 1000 flashes, about 4.46% of I was 
converted into III. Suppose that after one flash, a constant fraction x of the precursor 
is converted, then after 1000 flashes, the fraction converted is (l — (1 — a;) 1000 ), and 
equating this with the experimentally determined fraction of 0.0446, yields a conversion 
per flash x of 4.56 x 10 _3 %. Therefore, to correct for the effect of depletion, the para- 
hydrogen enhancement must be divided by, 

1000 x 

-, r « 1.023. 3.37 

(1 - (1 -x) 1000 ) V ; 

3.4 Data processing 

An estimate of the purity of the two qubit (Werner singlet) state is, 

x S x F x V f x —J—, (3.38) 



T x 1.023 1 30864' 

where P is some kind of intensity of the enhanced para-hydrogen signal and T is some 
measure of intensity of the thermal signal acquired from the multi-scan, multi-flash 
control experiment. This section considers the data processing steps we employed in 
interpreting the spectra and determining the appropriate values of P and T for the 
purity calculation and state tomography. 

3.4.1 Peak integration 

The enhanced and control spectra were processed |1911 1192j using home-written soft- 
ware 9 and analyzed by integration: measuring areas between the absorptive Lorentzian 
peaks and a flat, zero mean baseline. The integrals are a more reliable measure of sig- 
nal intensity than line heights due to two reasons. First, line heights critically depend 

Nuclear Overhauser Effect during inter-scan delays. 

9 The data processing was done on a Sun Blade 1000 machine running Solaris OS. Spectra were 

acquired on a Bruker spectrometer at the University of York, and the data was then converted and 

adapted for use on the Blade machine. I wrote the programmes and scripts for data conversion, 

manipulation, J matching, doubling and integration in C and in AWK 193 . Other processing tasks 

such as zero-filling, Fourier and inverse Fourier transformations and plotting/display were performed 

using software written by my supervisor. Finally, offset and baseline correction were done with 

Mathematica and the Grace [194) data visualization package. 
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on the frequency resolution |191| I192| of the spectrum, which is quite poor for our 
raw spectra; poorly digitized spectra often resulting in underestimates of the height. 
Second, the total integral under a band-unlimited frequency domain spectrum is equal 
to the first point in the time domain FID (besides a normalization factor) P95|, and 
is therefore independent of RF inhomogeneities and T 2 relaxation effects. 

However, prior to integration, it is important to deal with some spectral artifacts. A 
baseline offset or a sloping baseline introduces errors; the effect is illustrated in Figure 
I3.1U1 I employed second order curve fitting (baseline correction) prior to integration. 
Other important factors to consider are the need for properly phased, absorptive mode 
spectra and the number and extent of data points |191j . over which the peaks are 
integrated. The integrals also fluctuate due to noise and standard errors on the integral 
values were measured by calculating the variance over regions of baseline noise of 
identical spectral widths. 

3.4.2 J processing 

In the para-hydrogen experiment, the integrals of the anti-phase multiplet partially 
cancel due to overlap. Any peak integrals we directly measure, even after removing 
all the artifacts discussed above, will therefore be an underestimate. There is however, 
one way of reducing the overlap by "dragging" apart the peaks; the method is called J- 
doubling [1961 1197] . In summary, the data processing steps involved in accurate peak 
integration comprised the following steps. (I have described the reduction in peak 
integrals due to overlap and J-doubling in considerable detail in Appendix ?? and here 
I give only the major results.) 

1. We excised a small section from the overall spectrum that contained the dihydride 
peaks. 

2. The excised spectrum was inverse- Fourier transformed to yield a time domain 
signal and zero-filled |1911 1192j to increase the apparent frequency resolution. 

3. We implemented a computer algorithm scanning trial J' values, looking for a 
global minimum in the integrals of absolute J'-modulated spectra. This step 
results in an estimate for the true splitting, J. A plot of integrals of the J'- 
modulated spectra is shown in Figure 13.111 
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Figure 3.10: Effect of a baseline zero-offset on the peak integrals. The spectrum in (a) 
comprises an ideal anti-phase peak (maximum and minimum are between rs ±0.05) 
and (b) is the corresponding integral (obtained by summing up all previous points, 
starting from the left edge of the spectrum). The integral is uniformly zero in the 
baseline region, reaches a maximum « 0.5 in the region between the peaks and drops 
back to zero after the second peak. In (c), the spectrum (a) is raised upwards by 
a small amount (~ 0.001), and (d) shows the corresponding integral, indicating that 
now, even the baseline region has a positive integral which additively interferes with 
the peak integral. In this case, the peak integrals will be an over-estimate and baseline 
correction becomes necessary. 
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Figure 3.11: Results of the computer programme searching for the global minimum in 
the integrals of absolute J' modulated spectra (Appendix ??). The trial values, J' are 
given in units of the Nyquist frequency. The global minimum of the integral plot is 
indicated by an arrow and corresponds to the matched J' = J = 0.0054 value. 

4. The matched J value is then used to double the original spectrum. We double 
four times, substantially reducing the effects of overlap but also degrading the 
signal-to- noise ratio. The successively doubled spectra are shown in Figure EH21 

5. Peak integrals are determined from the doubled spectra and rescaled, taking 
into account the number of doublings employed. These rescaled integrals are 
given in Table 13.31 As expected, after doubling a sufficient number of times, 
the overlap disappears and the integrals no longer destructively interfere. As 
a result the measured peak integrals increase and the effect is shown in Figure 
13.131 Another important consequence of the doubling is to reduce the artificial 
imbalance between the peaks within a doublet. We take the average between 
the statistically indistinguishable peaks from the m = 4 spectra, resulting in 
accurate I and S peak integrals, which in this experiment are 8.0310 ± 0.1936 
and 7.4328 ± 0.1936 (calculated from the entries in the last row in Table EOj) . 



100 



(a) 



(b) 



(c) 



(d) 



(e) 



Figure 3.12: Effect of J-doubling. Only the left (S*) peak is shown. Part (a) shows 
the original spectrum (m = 0), and (b)-(e) show the successively doubled spectra 
(m = 1,2,3,4). The spectra are drawn to the same vertical scale and the signal 
heights are seen to drop by roughly half at each stage. 
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S, 



SD 



5.9381 6.8250 6.3643 6.8375 0.075 

1 6.6281 7.3230 7.0029 7.3095 0.128 

2 7.1190 7.2129 7.5793 7.2402 0.136 

3 7.2967 7.0603 8.0698 7.3507 0.090 

4 7.6220 7.2544 8.0594 8.0027 0.194 

Table 3.3: Integrals calculated from the spectra, where m is the number of J-doublings 
employed. The subscripts "1" and "r" refer to the left and right lines in the / and S 
spin peaks. 
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Figure 3.13: Peak integrals as a function of the number of J-doubling steps, m. The 
plot (a) shows the integrals for the S spin and (b) shows the integrals for the / spin. 
The integrals, in general, increase with increasing m and for m = 4, the peak integrals 
within a doublet become statistically indistinguishable. 
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3.4.3 State tomography results 

The I and S integrals of the processed spectra are then normalized by multiplying by 
the factor (S x F x V f )/(T x 1.023 x 30864), where T is the "raw" average integral 
under each of the four peaks, determined after baseline correction but without any 
J-processing. In our case, we measured a T of 303.12 ± 4.44. The other experimental 
parameters were: S = 3072, F = 1000 and Vf = 12.5/34 = 0.368. This calculation 
gives us the p and q coefficients in the state representation (J3.3)) . therefore achieving 
state tomography. However, we prefer to express the state in the form, (|3.2|) and use 
the equations ()3.6J) - (j3.7|) to determine the singlet and triplet fractions (a, b, c) from p 
and q. For example, in the particular experiment I am discussing, the / spin normalized 
integral corresponds to q — —8.0310 ± 0.1936 and the 5* spin integral corresponds to 
p = —7.4328 ± 0.1936. For a pure So, we expect p = q = — 1, however in our case, the 
multiplets are slightly imbalanced, resulting in unequal fractions of the triplets. Using 
(|3.6|) - (|3.7|) . the exact state characterization is, 

a = 0.9371, (3.39) 
b= 0.0448, and (3.40) 
c = 0.00905, (3.41) 

not identical but very similar to a Werner singlet with balanced triplets. The resulting 
density matrix is portrayed in Figure 13.141 and has a concurrence and entanglement of 
formation of 0.874 and 0.822 respectively. 

If we assume that our quantum state is a Werner singlet, (1 — 5)14/ 4 + eSq, with a 
singlet fraction equal to the a calculated above, then the polarization of the state is, 

e = 0.916 ±0.2, (3.42) 

which I call the effective purity and is consistent with a spin temperature of 6.4 mK 
or an effective magnetic field of 0.45 MT at room temperature. 

3.5 CW experiments 

This section discusses our isotropic PASADENA experiment utilizing the CW laser. 
The spectra from extended hydrogenation periods of 60 ms and 300 ms, in the presence 
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(a) (b) 

Figure 3.14: Results from quantum state tomography, showing the density matrices 
from (a) the state generated by the pulsed laser, compared with (b) a perfect singlet. 
Despite small non-zero amounts of T±\ in (a), the overall density matrices are very 
similar. The elements in the bar chart represent magnitudes of the elements in the 
density matrix written in the conventional {|a/3), \a/3), \/3a), |/3/3)} order. For example, 
the top left element in the graph represents the magnitude of the element (aa\p\aa) 
and so on. All elements in the matrix are real. 

of the MLEV-16 isotropic mixing sequence are shown in Figure 13.151 These spectra 
show that the doubly anti-phase character of the spectrum is maintained, indicating 
that the mixing sequence is very good at preserving the form of the singlet. However, 
there is one problem: naively one would expect the 300 ms spectrum to be five times 
bigger than the 60 ms spectrum, but it is not quite so large. We attribute this discrep- 
ancy to the relaxation of the highly polarized state, which can cause pronounced signal 
loss for extended durations of hydrogenation. A rough calculation gives a value of 
about 0.6 s for the relaxation constant. It is important to note that this value does not 
correspond to the single quantum T2 relaxation normally quoted in NMR, but rather it 
corresponds to a complicated decay of an isotropic state towards the maximally mixed 
state in the presence of a mixing sequence; this is sometimes described as "spin-lattice 
relaxation in the rotating frame", Tip |187j . Quantum state tomography results are 
outlined in Table 13.41 and the density matrices for the CW experiment are depicted in 
Figure 13.161 
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Figure 3.15: Spectra from the isotropic PASADENA experiment, employing MLEV-16 
isotropic mixing. The spectrum in (a) is acquired after 60 ms of laser irradiation and 
(b) shows the spectrum obtained after 300 ms of irradiation. Naively one would expect 
(b) to be five times bigger as compared to (a), however the signal intensities are smaller 
than expected. For comparison, (c) is a plot of the spectrum (a) multiplied by five. 



60 ms 300 ms 



a 


0.9159 


0.8191 


b 


0.0565 


0.0687 


c 


0.0138 


0.05615 


e 


0.8879 


0.7588 



Table 3.4: Tomography results from isotropic PASADENA experiments, involving CW 
irradiation for 60 and 300 ms; e is the effective purity. 
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Figure 3.16: Tomography results from the CW data, showing the density matrices after 
(a) 60 ms and (b) 300 ms irradiation. Build-up of significant T±\ error terms is evident 
from (b). Interpretation and labelling of this figure are identical to those in Figure 

EH 



106 



Chapter 4 

Implementation of quantum 
computation with the pure state 
quantum computer 

Chapter 01 described the experimental preparation of a two qubit pure quantum state — 
the qubits are the two X H nuclei in the organometallic compound Ru(H) 2 (C0) 2 (dppe), 
and the nuclei inherit the nuclear spin singlet state, derived from pure para-hydrogen. 
At this juncture, the natural question to ask is, having prepared a pure state, can we 
go one step further and actually use this as an initial state for a quantum information 
processing task, such as for implementing a quantum algorithm? The answer is yes 
and this chapter describes our experiments achieving this. 

The extremely high purity of the system, quite naturally, obviates the need for 
assembling pseudopure states, a process requiring non-unitary operations and generally 
lengthy or complicated pulse sequences. Our experiments [SHI I154j constitute the 
first implementation of quantum algorithms in liquid state NMR using pure initial 
states. As a result, they do not prepare pseudopure states, instead they initialize the 
system directly in an essentially pure quantum state (also lying above the entanglement 
threshold). In related work, Hiibler et al. implemented |152j the Deutsch algorithm on a 
para-hydrogen derived two qubit system, although with much smaller purities, e ~ 0.1. 
Our implementations extend their results, to almost perfectly pure quantum states, 
with e very close to 1. Furthermore, in addition to achieving high purity and the 
possibility for distilling entanglement [135J, our experiments are unique in the sense 
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System NMR parameters 
Ru(H) 2 (CO) 2 (dpae) 



resonance frequencies 



—7.61 and —7.22 ppm 

160 Hz (on the 400 MHz spectrometer) 



S 




4.8 Hz 



0.67 s 



1.9 s 



Table 4.1: Important NMR parameters of Ru(H) 2 (CO) 2 (dpae) used in the implemen- 
tation of Grover's search. 

that they prepare pure states on demand; our laser flash acts as a convenient switch, 
enabling us to generate the initial state, as and when desired. 

In this chapter, I describe the implementation of two quantum algorithms using pure 
initial states. They are solving Deutsch's problem of distinguishing between constant 
and balanced functions and performing Grover's search on an unstructured database of 
four quantum states. For the two algorithms, we employed two different molecular sys- 
tems: for Deutsch's algorithm, we used Ru(H) 2 (C0)2(dppe) — a chemical system with 
which we are quite familiar by now — and for Grover's search, we used its closely related 
arsenic analogue, Ru(H) 2 (CO) 2 (dpae), where dpae is l,2-bis(diphenylarsino)ethane. 
The important NMR parameters of the latter compound, whose structure is also shown 
in Figure 12.31 are summarized in Table 14.11 

Comparing the NMR parameters for Ru(H) 2 (CO) 2 (dpae) with Ru(H) 2 (C0) 2 (dppe) 
(see Table l?T2"|) . two results are obvious: first, the Ru(H) 2 (CO) 2 (dpae) has a slightly 
longer T 2 relaxation time, resulting in better decoherence properties and second, the 
frequency separation is about three times smaller in the "dpae" than in the "dppe" and 
this puts more stringent requirements on spin selective pulse sequences in the "dpae" 
system. These properties are important in evaluating the suitability of a molecu- 
lar system for NMR quantum computing, and are captured by the ratios (5/ J) and 
(J/(l/T 2 )): a system is suited for quantum information processing, if the condition 
6 ^> J ^> 1/T 2 is satisfied. These characteristic ratios for the two molecules are com- 
pared in Table 14.21 indicating that the "dpae" has slightly better decoherence and 
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Pfl T£l TY1 C*V 


"Hnnp" 


"Hnflp" 


S/J 


106.96 


33.33 


J/(l/T 2 ) 


2.67 


3.22 



Table 4.2: Frequency separation, splitting and decoherence parameters for the "dppe" 
and "dape" compounds. 

worse frequency selectivity properties. 

However, there is one clear cut advantage of using the Ru(H) 2 (CO) 2 (dpae) system: 
the initial state can be prepared with a purity indistinguishable from one; our mea- 
surements indicate a purity of e ~ 1.06 ± 0.2. Thus, we can describe our initial states 
as being pure, without further qualification. 

This chapter is organized as follows. Section 14.11 describes the Deutsch algorithm 
and Section l4~2l takes up Grover's search with the pure state quantum computer. Ex- 
perimental imperfections are visible in our spectra, and these can be roughly explained 
using a simple model based on representing quantum operations in terms of a family 
of Kraus operators. Appendix ?? exemplifies the use of this model, in context of the 
Deutsch spectra. The following two sections, more or less, reproduce our published 
and submitted works, and |154j . and so include brief discussions of many topics I 
have already discussed in previous chapters. 

4.1 Deutsch's algorithm with pure quantum states 

A published version of this Section can be found on the e-print server: 
http:// www.arxiv.ore; /quant-ph/0406044 



4.2 Grover's quantum search with pure quantum 
states 

A published version of this Section can be found on the e-print server: 
http:/ /www. arxiv.org/quant-ph/0407091 
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Chapter 5 

Practical implementation of twirl 
operations 

Entanglement is an important resource in the quantum setting and has been introduced 
in Chapter [T] Communication protocols, such as quantum teleportation |109t IllUj . 
work optimally, when they use maximally entangled pure states. For two qubits, these 
states are well known and are simply the four Bell states, given as = ( 1 00) ± 
|ll))/v^2 and |^ ± ) = (|01) ± |10))/\/2. In the standard teleportation experiment, 
Alice and Bob share a pair of maximally entangled pure states and the pre-existing 
entanglement, in fact, constitutes a noiseless quantum channel. In real life situations 
however, the channels are never perfect. The pure states are subject to decoherence and 
degrade into mixed states, entailing a decrease in both the purity and the entanglement. 
Entanglement concentration or purification protocols |202[ I203| then pick up these 
degenerated states and recover both the lost purity and entanglement, accompanying 
a decrease in the ensemble size |2U4j . The first step in most of these protocols is the 
"twirl" operation which converts an arbitrary mixed state p', not maximally mixed, 
into the Werner state p e , (jl.69j) . A useful fidelity measure (2021 for such states is their 
most significant overlap with the maximally entangled state, 

F{p) = max(^|p|^), (5.1) 

where the maximum is taken over all four Bell states. It is well known that entan- 
glement can be distilled from the mixed p whenever the state fidelity F > 1/2 [202J. 
Without loss of generality, we can also assume that the dominant term in p is the 
singlet state the fidelity F then becomes the overlap of the state with the singlet, 
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which I have referred to as the singlet fraction, (|1.84j) . An equivalent parametrization 
of p £ can thus be made in terms of F, 

p w = F\ip~)(ip~\ 

1_ F (5-2) 

+ -^-(^ + ><^ + i + + irxri), 

the correspondence between the polarization e of the Werner singlet and its singlet 
fraction being given as 

F(Ps) = (5.3) 

It can be checked that for the maximally mixed state e = and F — 1/4; whereas for 
the pure singlet e = 1 and F = 1. 

Consider a mixed state p' which satisfies two properties: first, it is not maximally 
mixed and second, it has a non-zero singlet fraction F. The twirl operation would take 
such a state to p £ , while conserving the singlet fraction, 

p' — > p £ , such that (5.4) 
F(p £ ) = F(p'). (5.5) 

The restraints on p' simply arise from the fact that maximally mixed states remain 
unchanged under all quantum operations, and that p' must have a non-zero singlet 
fraction because the twirl sequence does not create new singlet, rather it only preserves 
the singlet, while averaging out all other terms, converting them effectively to the 
maximally mixed state. In other words, the twirl sequence "Wernerizes" the input 
state, by letting the singlet pass unscathed, while scrambling other terms into the 
maximally mixed state. 

We have already come across the partial twirl in Section I3~T1 in context of preparing 
a state for one-shot tomography. This chapter discusses the full twirl: the different 
ways of practically implementing the operation, especially on an ensemble quantum 
computer. Section ?? introduces the implementation. Twirling is a kind of state 
averaging and Section ?? discusses this concept for a single qubit while Section ?? 
extends these ideas to the two qubit case. Finally, a step-by-step NMR implementation 
is discussed in Section ?? and then I conclude in Section ??. This material is closely 
based on our submitted manuscript, [20 5j. 

An updated version of this chapter can be found on the e-print server: 
http:// www.arxiv.org /quant-ph/0 409142 
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Chapter 6 
Purity sharing 



Preceding chapters dealt with different aspects of our para-hydrogen derived two qubit 
NMR quantum computer. Scaling up the prototype computer to larger-sized systems is 
important if we wish to address more challenging and practically useful mathematical 
problems. In this spirit, one can come up with different ideas for extending the size 
of our quantum computer. Some of these ideas seem promising but may prove quite 
difficult to implement in practice. 

1. For example, we can envisage a molecule (like an unsaturated organic polymer) 
which has multiple sites where H2 molecules can add forming a 2M qubit quantum 
computer, M being the number of dihydrogen molecules. 

2. Another scheme can be visualized, involving the transfer of purity from the two 
pure hydrogen qubits onto other weakly polarized qubits. After the polarization 
transfer, the hydrogens are allowed to detach from the molecule and their place is 
taken up by fresh, pure qubits and the transfer step is repeated. In this manner, 
purity can be "pumped" into the molecule. This scheme to some extent resembles 
the method of algorithmic cooling [96J already mentioned in Section 11.2.41 

3. Both of the above ideas involve bringing purity into the molecule from outside. 
We can also think of another scenario: suppose we are given a molecule with 2M 
pure or nearly pure qubits, derived from M para-hydrogen molecules. Instead 
of "importing" more purity into the molecule, we can consider the question of 
sharing out the purity on k > 2M qubits. The question that I investigate in 
this chapter is whether this sharing can be carried out, while remaining above the 
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entanglement threshold. This question is related to purity compression strategies 
outlined by Schulman and Vazirani [M], but while the latter schemes concentrate 
purity onto a smaller number of spins, what I consider here is "diluting" or 
sharing out this purity onto a larger quantum subspace. 

An updated version of this chapter can be found on the e-print server: 
http:// www.arxiv.ore; /quant-ph/0509036 
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